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THE SPECIFICATION OF SUGARSCAPE 


JOSEPH KEHOE 


Abstract. Sugarscape is a well known and influential Agent Based Social Sim¬ 
ulation (ABSS). Various parts of Sugarscape are supplied as examples in almost 
all Agent Based Model (ABM) toolkits. It has been used for demonstrating 
the applicability of different approaches to ABM. However a lack of agreement 
on the precise definition of the rules within Sugarscape has curtailed its useful¬ 
ness. We provide a formal specification of Sugarscape using the Z specification 
language. This demonstrates the ability of formal specification to capture the 
definition of an ABM in a precise manner. It shows that formal specifications 
could be used as an approach to tackle the replication problem in the field of 
ABM. It also provides the first clear interpretation of Sugarscape identifying 
areas where information is missing and/or ambiguous. This enables researchers 
to make proper comparisons between different implementations of this model. 


1. Introduction 

1.1. Overview. First we give an short informal snmmary of Sngarscape. We 
follow with a brief introduction to formal specihcation. We then specify the single 
resource simulation. Following the standard Z patterns of development we list 
the basic types and constants hrst, followed by the specihcation of the basic state 
attributes and invariant properties. Then the rules are presented in order. 

After this we present the extended specihcation, that is the specihcation with 
two resources. Here we highlight the diherences between the single resource and 
two resource specihcations by presenting the specihcations in bold face whenever 
the specihcation changes from the original. 

2. Sugarscape 

2.1. Agent Based social Simulations. Sugarscape was the hrst large scale 
Agent Based Social Simulation (ABSS). It was developed by Epstein and Axtell 
and presented in their book Growing Artificial Societies [Epstein and Axtell, 1996]. 
The release of this simulation is considered an important event in the emerging 
held of Agent Based Social Simulation. 

The Sugarscape ABSS was used to investigate how individual behaviour can 
inhuence and cause diherent social dynamics within large populations. It has been 
used to show how, for example, inheritance of wealth ahects resource distribution 


Date: 9/03/2015. 


1 



2 


JOSEPH KEHOE 


in populations and how disease can spread through a population. It remains influ¬ 
ential today and every major simulation toolkit (Swarm, Repast, Mason and Net- 
Logo) [Railsback et al., 2006, Berryman, 2008, Inchiosa and Parker, 2002] comes 
with a partial implementation of Sugarscape that demonstrates that toolkit’s ap¬ 
proach to simulation. Since Sugarscape first appeared ABSSs have been applied 
to fields as diverse as Anthropology[Campillo et ah, 2012], Biomedical Science, 
Ecology, Social Science [Axtell and Axtell, 2000], Epidemic modelling and Market 
Analysis[Macal and North, 2009, Troitzsch, 2009, Gilbert, 2004]. 

ABSS’s employs a bottom-up approach to modelling populations. Instead of pre- 
computing the overall population behaviour, as done in equation based models, 
individual agents and their local interactions within the population are modelled. 
The behaviour of the overall population is left to emerge from these local inter¬ 
actions. This approach allows us to address failings in the top-down approach 
and demonstrates the causal factors behind the emergence of group dynamics. In 
cases where we do not know what the overall behaviour will be or where we are 
trying to And out the causes of this behaviour, bottom-up based ABSSs are the 
only possible approach. 

2.2. Issues with Sugarscape. Currently, social science simulations are start¬ 
ing to embrace concurrency in an effort to allow for bigger, more complete and 
faster implementations of ABMs. Different concurrency researchers have used the 
Sugarscape model as a testbed for benchmarking different approaches to parallelis¬ 
ing ABMs [Lysenko and D’Souza, 2008, Perumalla, 2006, Richmond et al., 2010]. 
However although the rules of Sugarscape have been defined there is no general 
agreement on their exact meaning. These difficulties hamper the ability of re¬ 
searchers both to properly compare their approaches, provide complete implemen¬ 
tations of Sugarscape or replicate their results. 

Most of the rules require some form of conflict resolution. We have specified 
the rules in a manner consistent with the original intention (agents acting con¬ 
currently) but independent of any particular approach to how this concurrency 
is implemented. That is, we have refrained from imposing any specific conflict 
resolution rules. 

By formalising Sugarscape and providing a single precisely defined reference for 
the rules we can produce a standard definition of Sugarscape. Compliance with 
this single reference will allow proper comparisons to be made between different 
approaches. It also leaves it open to the implementer to decide what approach to 
conflict resolution they wish to take. We detected ambiguities present in the cur¬ 
rent rule definitions, provided precise interpretations, where possible, and flagged 
irresolvable problems where not. 

We made the decision to restrict the initial specification to one pollution type 
and one resource type in an effort to guarantee clarity. While the rules were 
designed so that they could be extended to arbitrary numbers of resources and 
pollutants, explicitly specifying for an arbitrary number of resources and pollutants 
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would make the specification even more difficult to understand and thus more likely 
to either contain or cause mistakes. 

Once we had a specification for the single resource scenario we extended the 
specification to a two resource situation. This allowed us to specify the final rule, 
Trade, as that rule requires two resources to function. 

This allowed for: 

(1) A simpler and easier to understand specification of the rules that use only 
one resource (trading clarity against completeness); 

(2) A complete (but separate) specification for simulations that use two re¬ 
sources. 

We do not provide specifications for multiple pollutants as multiple pollutants were 
never actually implemented in Sugarscape^. 

Similarly we did not provide a specification for more than two resources as we 
deem the benefits of doing so counterbalanced by both the complexity of the result¬ 
ing specification and the lack of any requirement to use such a complex simulation 
for benchmarking purposes. Sugarscape has only ever been implemented with two 
resources types, known respectively as sugar and spice. Anyone wishing to extend 
Sugarscape further can use the two resource specification for guidance. 

2.3. Synchronous and Asynchronous Updating. Originally the rules were 
stated with an explicit assumption that the underlying implementation would be 
sequential. Concurrency was simulated through randomisation of the order of each 
rule application on the individual agents, and models that follow this regime are 
termed asynchronous. 

All results reported here have been produced by running the model 
on a serial computer; therefore only one agent is “active” at any 
instant. In principle the model could be run on parallel hardware, 
permitting agents to move simultaneously (although M would have 
to be supplemented with a conflict resolution rule to handle cases in 
which two or more agents simultaneously decide to inhabit the same 
site). 

[Footnote 12, Chapter II] 

The alternative to asynchronous updating is synchronous updating. Synchro¬ 
nous updating assumes that all updates occur concurrently. While it is clear that 
the original authors have no objection to employing synchronous updating on sug¬ 
arscape it is well known that asynchronous and synchronous updating produce 
different results. What is not known is how divergent these results are in the 
case of complex ABMs such as Sugarscape or indeed how to apply synchronous 
updating to all the complex interaction types in Sugarscape. 


^We leave this open as an exercise for the reader. 
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In order to answer these questions we present initially a specification that as¬ 
sumes a synchronous updating regime, as this is the most novel approach. Follow¬ 
ing this we give the equivalent Asynchronous updating version for comparison. 

3. Single Resource Sugarscape 

Sugarscape is a discrete turn based simulation composed of a set of interacting 
agents that move across an environment. The environment, or simulation space, is 
modelled as a two dimensional M by M grid or matrix of discrete locations known 
as the lattice. This lattice is toroidal in nature, that is, it wraps around on all four 
edges. Every lattice location has a position denoted by its x and y coordinates. 
For any lattice location [i, j] there are four direct (von Neumann) neighbours (up, 
down, left and right) at positions [i, (j -t- 1)%M], [i, (j — 1)%M], [(i — l)%M,j] and 
[(i -|- l)%M,j]. We denote this set of von Neumann neighbours as Ni{i,j), and 
further use Nk{i,j) to denote the set of von Neumann neighbours where each is a 
maximum distance of k locations from the location [i,j]. 

Each location can hold a number of resource and pollution types. While there is 
no limit placed on how many resource or pollution types can exist in a Sugarscape 
simulation we are unaware of any Sugarscape derived simulation that uses more 
than two resource types and one pollution type. When there is only one resource 
type it is called sugar and if there are two then the second resource is known as 
spice. These amounts are measured as natural numbers (> 0). Each individual 
location has limits placed on the maximum amount of resources of any type it may 
carry at any one time. These limits are dehned at simulation startup and remain 
hxed during a simulation run. Agents consume the resources at their current 
location. Locations replenish their resources by some dehned amount during each 
time step. Each location can also hold at most one agent at a time. 

Agents reside at locations within the lattice but are mobile and can change 
location at most once per step. At a minimum each agent has the following 
attributes: 

Metabolism Rate (one per resource type): The rate at which an agents 
resource stores decrease during each simulation step. Different resource 
types have independent metabolism rates. Once an agent runs out of re¬ 
sources it dies (is removed from the simulation); 

Age: The number of steps that the agent has been present in the simulation; 

Maximum Age: The maximum number of steps that an agent is allowed to 
exist during the simulation run. Once an agent reaches its maximum age 
it is removed from the simulation; 

Resource Store (oue for each resource type): The amount of each re¬ 
source that the agent currently has; 

Visiou: How far in each of the cardinal directions that the agent can see. 
An agent can only interact with locations and agents that are in its neigh¬ 
bourhood N.^ision- To ensure locality all agent values for vision will be less 
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than some predefined maximum and this maximum will be much smaller 
than the lattice dimension size (M). 


In the more complex versions of Sugarscape agents can also have a “culture” 
identifier (identifying which tribe the agent belongs to), a set of outstanding loans 
of resources that the agent has given to (or received from) other agents, a set of 
diseases that the agent has contracted and, an immunity system that gives each 
agent immunity from certain diseases. 

A simulation run consist of a series of turns or steps during which certain rules 
are applied to each location and agent. Each rule is applied concurrently and 
instantaneously to each agent and/or location. The rules are generally fairly simple 
and the only information that an agent (or location) can use when deciding how 
to apply a rule is local information, that is an agent or location at position [i,j] 
can only access information from locations and/or agents that are within the set 
Nk, where k< vision (in most cases k=l). 

The rules for locations decide how resources are replenished and how pollution 
is created or spread. The rules for agents are more varied and determine agent 
movement and interaction. Agent interaction can range from spreading disease, 
trading, entering financial agreements and even combat. There are a large number 
of rules but not all rules need to be (or indeed can be) applied in the same simula¬ 
tion run. The rules are chosen based on what we wish to model. A simulation that 
wishes to see the effect of trading on wealth distribution would have no need for 
the combat or culture rules while one modelling disease transmission would only 
require the movement and disease transmission rules. 


3.1. Basic Types and Constants. First we identify the basic types and any 
required constants. Many are self explanatory or will become clear when their 
associated rules are specified. A simulation is defined by the values given to these 
constants and the combination of rules employed. 
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RATE :A 

INITIALSUGARMIN, INITIALSUGARMAX 
WINTERRATE, SEASONLENGTH : Ni 
PRODUGTION, GONSUMPTION : N 
GOMBATLIMIT : N 
IMMUNITY LENGTH : N 
INITIALPOPU RATION SIZE : N 
POLLUTION RATE : M 
OHILDAMT : N 


OULTUREOOUNT : Ni 
MAXVISION : Ni 

MINMETABOLISM, MAXMETABOLISM : N 

SUGARGROWTH : Mi 

MAX AGE, MIN AGE : Mi 

MAX SUGAR : Mi 

DURATION : Mi 


M ; M 


1 


1 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 
(7) 

(8a) 

( 86 ) 

M(9) 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 


OULTUREOOUNT mod 2 = 1 
MINMETABOLISM < MAXMETABOLISM 
MAX AGE < MIN AGE 
MAXVISION < M 

INITIALSUGARMIN < INITIALSUGARMAX 
INITIALPOPU RATION SIZE <M*M 


(1) The simulation space is represented by a two dimensional M by M matrix 
of locations. Each location in the simulation space is referenced by two 
indices representing its position in this matrix; 

(2) OULTUREOOUNT determines the size of the bit sequence used to rep¬ 
resent cultural allegiances. This is always equal to an odd number so that 
the number of I’s in the sequence is never equal to the number of O’s; 

(3) Agents can only “see” in the four cardinal directions, that is the locations to 
the north, south, east and west. Agents are endowed with a random vision 
strength that indicates how many locations the can “see” in each direction. 

This endowment is always less than MAXVISION and MAXVISION 
is always less than M; 

(4) Agents consume an amount of sugar (resources) during each turn. This 
sugar represents the amount of energy required to live. Each agent is en¬ 
dowed, on creation, with a random metabolism between MINMETABOLISM 
and MAXMETABOLISM] 

(5) Agents consume sugar (resources) from the location they occupy. Each 
location can renew its sugar at a rate determined by SUGARGROWTH. 
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After each turn up to a maximum of SUGARGROWTH units of sugar 
are added to each location (in accordance with the Growback rule); 

(6) MAX AGE and MIN AGE are, respectively, the maximum and minimum 
allowable lifespan for any agent; 

(7) MAX SUGAR is the maximum amount of sugar that any location can 
possibly hold. This is known as the carrying capacity of a location; 

(8) RATE and DURATION are used for determining the rate of interest 
charged for loans and the duration of a loan; 

(9) INITIALSUGARMIN and INITIALSUGARMAX are the lower and 
upper limits for initial endowment of sugar given to a newly created agent; 

(10) If seasons are enabled then two seasons, winter and summer are allowed 
with a duration of SEASON LENGTH turns (ticks) and a new separate 
lower seasonal grow back rate calculated using WINTERRATE (as de¬ 
termined by the SeasonalGrowback rule); 

(11) Pollution can occur at a rate determined by the production and consump¬ 
tion of resources determined by the PRODUGTION and GONSUMPTION 
constants respectively; 

(12) The combat rule posits the maximum reward GO MB AT LI MIT that can 
be given to an agent through killing another agent; 

(13) Immunity in agents is represented using a fixed size sequence of bits of 
length IMMUNITY LENGTH- 

(14) We have some predetermined initial population size INITIALPOPULATIONSIZE 
that is used to initialise the simulation; 

(15) POLLUTIONRATE determines the number of steps that elapse before 
pollution levels diffuse to their neighbours; 

(16) A certain amount of sugar reserves, OHILDAMT, are required for an 
agent to have children. 


[AGENT] (1) 

POSITION == 0 .. M - 1 X 0 .. M - 1 (2) 
SEX ;:= male \ female (3) 

BIT ;:= 0 I 1 (4) 

affiliation ::= red \ blue (5) 

boolean ::= true \ false 


(1) AGENT is used as a unique identifier for agents; 

(2) POSITION is also used to make specifying indices within the grid so as 
to make the schemas easier to read and more compact; 

(3) All agents have a sex attribute; 

(4) BITs are used to encode both culture preferences and diseases of agents; 
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(5) Every agent has a cultural affiliation of either belonging to the blue tribe 
or red tribe. 

Agents can, using the Mating rule, have offspring if they are fertile. Fertility is 
determined by the age of the agent, where fertility starts at some predefined age 
and ends at another. These boundaries are dehned for all agents. The numbers 
are set out by Epstein and Axtelland although there appears to be no special sig- 
nihcance attached to these numbers we will stick with the originals. Male fertility 
ends 10 turns later than female fertility. 

FEMALE FERTILITY START, FEMALEFERTILITYEND : N 
MALE FERTILITY ST ART, MALEFERTILITYEND : N 

12 < FEMALEFERTILITYSTART < 15 
40 < FEMALEFERTILITYEND < 50 
12 < MALE FERTILITY START < 15 
50 < MALEFERTILITYEND < 60 

MALEFERTILITYEND = FEMALEFERTILITYEND + 10 


The replacement rule requires a sugar allocation be given to new agents set 
between 5 and 25. Again there appears to be no special signihcance attached to 
these numbers. 

STARTSUGARMIN, ST ART SU G ARM AX : N 

STARTSUGARMIN = 5 
ST ART SUG ARM AX = 25 


3.2. The Sugarspace Lattice. The simulation space in Sugarscape consists of 
a finite discrete two-dimensional array of locations. Each location is identified its 
row and column value. Each location contains a number of resources. While only 
two resources are ever used it is clear that the intention of the original authors was 
that the simulation could be extended so that any number of different resources 
can be present. 

Similarly each location can contain a number of pollutant levels. In practice, 
although the rule is explicitly dehned for an arbitrary number of pollution types 
only one is ever used. Again, in line with actual Sugarscape usage and to make 
the specihcation more readable we assume only one pollution type. Pollution 
huxes are used in the rules to help calculate how pollution levels change over 
time. Although explicitly referenced in the Pollution rule these do not need to be 
explicitly modelled in the specihcation. 
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_ Lattice 


sugar : POSITION ^ N (1) 

maxSugar : POSITION ^ N (2) 

pollution : POSITION M (3) 


dom sugar = dommaxSugar = dompollution = POSITION (4) 

Vx : POSITION • sugar{x) < maxSugar{x) < MAXSUGAR (5) 


(1) sugar is a mapping that stores the amount of sugar stored at each position; 

(2) maxSugar is a mapping that records the maximum amount of sugar that 
can be stored (carried) in each position; 

(3) pollution records the amount of pollution at each location; 

(4) Every position has a sugar level, a maximum allowed sugar level (or carry¬ 
ing load) and a pollution level; 

(5) Every position’s sugar level is less than or equal to the maximum al¬ 
lowed amount for that position which is in turn less than or equal to the 
MAXSUGAR constant; 

We need to track the number of turns that have occurred in the simulation. 

Each turn consists of the application of all rules that form part of the simulation. 

Step _ 

step : N 


3.3. Agents. Every agent is situated on a location within the grid and each lo¬ 
cation is capable of containing only one agent at a time (putting an upper limit 
on the number of possible agents). Agents are mobile, that is they can move to 
a new location if a suitable unoccupied location is available. Movement is both 
discrete and instantaneous, it is possible for an agent to move to a new location 
instantly while skipping over all intermediate locations. The attributes that every 
agent has are: 

Vision: How far in the four cardinal directions that an agent can see; 

Age: Number of turns of the simulation that an agent has been alive; 

Maximum Age: Age at which an agent dies; 

Sex: Agents are either male or female; 

Sugar Level: The amount of sugar that an agent currently holds. There is 
no limit to how much sugar an agent can hold; 

Initial Sugar: The amount of sugar the agent was initialised with on cre¬ 
ation; 

Metabolism: The amount of energy, dehned by sugar (or resource) con¬ 
sumption, used during every turn of the simulation; 

Culture Tags: A sequence of bits that represents the culture of an agent; 

Children: For each agent we track its children (if any). To apply the In¬ 
heritance rule the full list of an agents children is required. 
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Loans: Under the credit rule agents are allowed lend and/or borrow sugar 
for set durations and interest rates so we need to track these loans. For 
each loan we need to know the lender, the borrower, the loan principal and 
the due date (represented as the step number); 

Diseases: Diseases are sequences of bits that can be passed between agents. 
An agent may carry more than one disease; 

Immunity: Each agent has an associated bit sequence that confers immu¬ 
nity against certain diseases. If the bit sequence representing a disease 
is a subsequence of an agents immunity bit sequence then that agent is 
considered immune to that disease. 
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_ Agents _ 

population : F AGENT 

position : AGENT ^ POSITION 

sex : AGENT SEX 

vision : AGENT Ni 

age ; AGENT ^ N 

maxAge : AGENT Ni 

metabolism : AGENT -H- N 

agentSugar : AGENT -H- N 

initialSugar : AGENT -H- N 

agentGulture : AGENT -H- seq i?/T 

children : AGENT AGENT 

loanBook : AGENT ^ [AGENT x (N,N)) 

agentimmunity : AGENT -^se(\BIT 

diseases : AGENT -t-)-Pseqi?/T 

population = 

dovaposition = domsex = dom msion 
= dom max Age = dom agentSugar = dom children 
= dom agentGulture = dom metabolism = dom age 


= dom agentimmunity = dom diseases 

(1) 

dom loanBook C population 

(2) 

dom[idmloanBook) Y population 

( 3 ) 

\/x,y: AGENT- d : seq BIT • 


x,y E population A x ^ y ^ 

( 4 ) 


[[age[x) < maxAge[x) A MIN AGE < maxAge[x) < MAXAGE 
A # agentGulture[x) = GULTUREGOUNT 
A 4k Immunity [x) = IMMUNITY LENGTH 
A ms*on(x) < MAXVISION 

A MINMETABOLISM < metaboUsm[x) < MAXMETABOLISM 
A position[x) = position[y) ^ x = y) 
d e ran diseases (x) ^ #d < IMUNITY LENGTH 


(1) Every existing agent has an associated age, sex, vision, etc. Note that the 
population holds only the currently existing agent IDs; 

(2) Only current members of the population can be lenders; 

(3) Only current members of the population can be borrowers 

(4) Every agent in the population is guaranteed to have a current age less than 
the maximum allowed age for that agent, a maximum age less than or 
equal to the global MAXAGE, a metabolism between the allowed limits 
and vision less than or equal to the maximum vision. The sequence of bits 
representing its culture tags is GULTUREGOUNT in size while those 
representing immunity is IMMUNITY LENGTH in size. All diseases 
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are represented by sequences of bits that are shorter than the immunity 
sequence. 

The entire simulation consists of locations, agents and a counter holding the tick 
count. We combine them all in the schema SugarScape. 


_ SugarScape 

Agents 

Lattice 

Step 


The initial state of the schema when the simulation begins must also be stated. 

_ InitialSugar Scape _ 

Sugar scape' 


step' = 0 (1) 

# population' = INITIALPOPULATIONSIZE (2) 

loanBook' = 0 (3) 

Va : AGENT • (4) 

a G population' ^ 


{age{a) = 0 A diseases' {a) = 0 A children' (a) = 0 

A INITIALSUGARMIN < agent Sugar'(a) < INITIALSUGARMAX) 
A initial Sugar'(a) = agentSugar' {a) 


(1) step is set to zero; 

(2) The population is set to some initial size; 

(3) There are no loans as yet; 

(4) Every agent in the starting population has an age of zero, no diseases or 
children and some initial sugar level within the agreed limits. The other 
attributes have random values restricted only by the invariants; 

3.4. Rules. There are a number of rules that can be employed in different com¬ 
binations to give different simulations. We will quote each rule as laid out in the 
appendix of [Epstein and Axtell, 1996] and follow, where necessary, with a more 
detailed explanation of the rule. In many cases the simple rule dehnitions are not 
complete. Extra information, embedded in the original text, has been extracted 
where necessary to help complete these rules. The majority of rule dehnitions 
assume only one resource (sugar) and it is these that are specihed in this section. 

The simulation is discrete with each time interval representing one complete set 
of rule applications. We use the step variable in the SugarScape schema to keep 
track of the current time interval number. 

Where there exist ambiguities in the rule dehnitions we will identify them and 
propose one or more possible interpretations consistent with what we believe to 
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be the authors intentions. Throughout the rule dehnitions constants such as a, jS 
are used but they have different meanings in each rule. For the sake of clarity we 
will give each constant a meaningful and globally unique name. 

3.5. Tracking Steps. While not dehned explicitly as a rule, we must ensure that 
we record the current step number. We increment the Step variable before every 
sequence of rule applications that compose a single turn of the simulation. 

There is an issue with metabolism in that every turn of the simulation requires 
that agents use up their sugar reserves at a rate determined by their metabolism. 
It is not explicitly stated when or where this sugar deduction occurs within the 
rules. It could be placed, for example, in the movement rule but it can also be 
placed, just as validly, within any rule that is guaranteed to be applied during 
every turn. Since there is no obvious reason why one is superior to the other, as 
long as it is consistently applied, we choose to place the metabolism deduction 
within the Tick schema. This new rule can be stated simply as follows: 

Tick: At the start of every time interval increase every agents age by one 
and decrease every agents sugar level by their metabolism rate. 


Tick _ 

AAgents 
A Step 

population' = population 
position' = position 
sex' = sex 
vision' = vision 
maxAge' = maxAge 
metabolism' = metabolism 
initialSugar' = initialSugar 
agentCulture' = agentCulture 
children' = children 
loanBook' = loanBook 
a gent Immunity' = agentimmunity 
diseases' = diseases 

step' = step +1 (1) 

V X : AGENT • X E population 

age'{x) = age{x) + 1 (2) 

A agent Sugar'{x) = agentSugar{x) — metabolism{x) (3) 


(1) Add one to the step count; 

(2) Increase everyone’s age by one; 

(3) Decrease everyone’s agentSugar by their metabolism. 
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3.6. Sugarscape Growbacka. 

Sugarscape Growback^: At each Lattice position, sugar grows back at a 
rate of a units per time interval up to the capacity at that position. 

Growback determines the rate at which location resources are replenished. The 
integer constant a indicates the amount by which resources grow during a single 
step or time interval. If a = cxo then each resource returns to its maximum value 
during each turn, i.e. it is instantly fully replenished after each step. The rule 
only refers to a single resource, sugar, but the book explicitly defines one other 
resource spice and it is clear that generalisations allowing an arbitrary number of 
resource types to be held at each Lattice position are acceptable. 

Since we are dealing only with one resource, sugar, we only need to dehne a for 
this resource . The constant SUGARGROWTH represents alpha in this rule and 
we use this to update the sugar level of each position. 

Since the maximum carrying level of each resource cannot be exceeded we will 
set the resource levels to its maximum value if application of the replenishment 
rate would result in a value greater than this maximum. With these dehnitions we 
can express the Growback rule in a simple manner. The last line in the schema 
(see below) does the work of updating the resource levels of every location. 

_ Growback _ 

ALattice 

pollution' = pollution 
maxSugar' = maxSugar 
sugar' = {x : POSITION • 

X niin{{sugar{x) + SUGARGROWTH, niaxSugar{x)})} (1) 

(1) The new sugar levels are calculated using a simple formula to either, 
the maximum possible level for that location or the old level plus the 
SUGARGROWTH whichever is the smaller. 

3.7. Seasonal Growback Sa^g^^. 

Seasonal Growback Initially it is summer in the top, half of the 

Sugarscape and winter in the bottom half. Then every 7 time periods the 
seasons flip - in the region where it was summer it becomes winter and vice 
versa. For each site, if the season is summer then sugar grows back at a 
rate of a units per time interval; if the season is winter then the grow back 
rate is a units per fl time intervals. 

Seasonal growback is an alternative to the previous grow back rule. Which rule 
is chosen will depend on what the simulation is trying to demonstrate. Seasonal 
grow back allow us to introduce seasonal factors into the original Growback rule. 
There are two seasons (representing summer and winter) and each lasts 7 turns 
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before switching. We rename 7 to SEASON LENGTH, a is the summer season 
SUGARGROWTH rate and a/(3 is the winter season rate. We use the existing 
SUGARGROWTH to hold the summer rate and introduce WINTERRATE as 

(3. 

Determining what season it is during a turn is fairly trivial. When seasonLength 
divides into the Step variable evenly it is summer in the top half and winter in the 
bottom half (and vice versa). 

_ SeasonalGrowback _ 

ALattice 

EStep 


pollution' = pollution 
maxSugar' = maxSugar 
Vx : POSITION* * 

{step div SEASONLENGTH) mod 2 = 0^ sugar' = (1) 

{x : POSITION I first{x) < M div 2 • 

X i-A min{{sugar{x) + SUGARGROWTH, maxSugar{x)})} (la) 

U 

{x : POSITION I first{x) > M div 2 • (16) 

X eA min{{sugar{x) + SUGARGROWTH div WINTERRATE, maxSugar{x)})} 
{step div SEASONLENGTH) mod 2^0^ sugar' = (2) 

{x : POSITION I first{x) < M div 2 • 

X ^ min{{sugar{x) + SUGARGROWTH div WINTERRATE, niaxSugar{x)])]{2a) 
U 

{x : POSITION-, y:n\ first{x) > M div 2 • 

X HA min{{sugar{x) + SUGARGROWTH, maxSugar{x)})} (26) 


(1) If the season is summer then: 

a) Top half of grid is updated as normal; 

b) Bottom half is updated at winter rate. 

(2) Otherwise if it is winter: 

a) Top half of grid is updated at winter rate; 

b) Bottom half is updated as normal. 

3.8. Movement - M. 

Movement - M: 

• Look out as far as vision permits in each of the four lattice directions, 
north, south, east and west; 

• Considering only unoccupied lattice positions, hud the nearest position 
producing maximum welfare; 

• Move to the new position 

• Collect all resources at that location 






16 


JOSEPH KEHOE 


The previous rules affected only the locations but the remaining rules affect 
agents as well as locations. The Movement rule determines how agents select 
their next location. There are a number of different versions of this rule. We will 
specify the simplest rule hrst as it is the only movement rule explicitly dehned in 
the appendix but we will also specify the other movement rules dehned in the text. 
We add a subscript to the rule title (Mbasic) to distinguish between the different 
movement rule specihcations. 

Not explicitly stated within the rule but stated as a footnote to the rule is the 
restriction that the order in which the lattice directions are searched should be 
random. This comes into play when two or more available sites exist with the 
same welfare score. 

This rule does not guarantee that an agent will move to the best location. 
To see why this is the case consider what happens if two agents both try to 
move to the same location. Only one can succeed and the other will have to 
move to a less advantageous location. How we decide which agent succeeds is 
not dehned. We assume that either a conhict resolution or conhict avoidance rule 
is available to make this decision but it is not stated what this rule should be. 
The original implementation is sequential with agents assumed to be moving in 
a random order thus enforcing collision avoidance. No guidance is provided for 
concurrent implementations. 

To help make the specihcation clear we dehne some simple helper functions. 
The distance between two positions is only dehned for positions that are directly 
horizontal or vertical to each other. This function must take into account the 
torus-like (wrap around) structure of the simulation. 


distance : POSITION 
X POSITION 


Vxl, x2, yl, y2 -.N • 
distancel(xl, yl), (xl,y2)) = (1) 

min{{\ y2 - yl |, M- | y2 - yl |}) 
distance{{xl,yl), {x2,y2)) = ( 2 ) 

min{{\ xl — x2 |, M— | xl — a;2 |}) 
distance{{xl,yl), {x2,y2)) = oo 

xl ^ x2 Ayl ^ y2 (3) 


(1) If two agents are vertically aligned we calculate distance based on the 
horizontal distance; 

(2) If two agents are horizontally aligned we calculate distance based on the 
vertical distance; 

(3) Otherwise the distance is dehned as inhnity. 
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We use this to define the adjacent function that lets us know if two agents are 
directly beside each other. 


adjacent : POSITION 
X POSIT ION 
-^boolean 

ya,b: POSITION • 
adjacent{a, b) ^ distance{a, b) = 1 


visibleAgents takes an agent, a function mapping agents to positions and the 
vision range of the agent and returns the set of agents that are within that agent’s 
neighbourhood. 


visibleAgents ; AGENT 
X {AGENT POSITION) 
xN 

AGENT 

y agent : AGENT] pos : AGENT POSITION] range : M • 
visibleAgents{agent,pos, range) = 

{ag : AGENT \ ag G dompos A 1 < distance{pos{ag),pos{agent)) < range} 
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_ M ovement basic _ 

ASugarScape 

step’ = step 

population' = population 

maxSugar' = maxSugar 

pollution' = pollution 

sex' = sex 

vision' = vision 

age' = age 

maxAge' = max Age 

agentCulture' = agentCulture 

loanBook' = loanBook 

diseases' = diseases 

agent Immunity' = agent Immunity 

children' = children 

metabolism' = metabolism 

initial Sugar' = initial Sugar 

Va : AGENT] I : POSITION • 


a G population' ^ (1) 

distance{position'{a),position{a)) < vision{a) 

{distance{position{a),l) < vision{a) A (/ ^ run position')) ^ (2) 

sugaril) < sugar {position' {a)) ( 2 a) 

A {distance{l,position{a)) < distance{position'{a),position{a))) (26) 
^ sugar{l) < sugar {position' {a)) 

agentSugar' = {Va : AGENT \ a G population' • 

a I—)■ agentSugar{a) + sugar {position' {a))} (3) 

sugar' = sugar © (V/ ; POSITION \ I G lanposition' 0} (4) 


After the rule is applied the following will be the case for every agent: 

(1) They will be located within one of the locations in their original neighbour¬ 
hood (possibly the same position as before); 

(2) After every agent has moved: 

a) There will exist no remaining available locations from the original neigh¬ 
bourhood of an agent that would have given a better welfare score than 
the location that agent now inhabits (we picked the maximum reward); 

b) If there was more than one location with maximum reward then the 
agent moved to the closest location. 

(3) Agent sugar levels increase because they consume all the sugar at their new 
location (even if the new location is the same as their old location); 
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(4) Location sugar levels are set to zero everywhere there is an agent present. 


The specihcation states what is true after the application of the rule but not how 
we achieve that state. In any implementation some conflict resolution strategy will 
be needed but in this specihcation we remain agnostic as to what it should be. 

The rule is well stated but requires that we precisely dehne welfare. For a single 
resource simulation welfare is precisely equal to the amount of sugar available at 
a location. We will dehne welfare for multiple resource simulations later. 


3.9. Pollution Formation Pn,x- 


Pollution Formation Th,/?: When sugar quantity s is gathered from the 
Sugarscape, an amount of production pollution is gathered in quantity as. 
When sugar amount m is consumed (metabolised), consumption pollution 
is generated according to f3m. The total pollution on a site at time t, p^, 
is the sum of the pollution present at the previous time, plus the pollution 
resulting from production and consumption activities, that is, p* = + 

as + (3m. 


This single resource pollution rule is easiest to understand and the most common 
form of the pollution rule. When pollution is incorporated into the Sugarscape the 
movement rule is changed so that the welfare of a location is now dehned using 
the sugar to pollution ratio - the greater the ratio the greater the welfare. This 
ratio is dehned as sugar/ (1 + pollution) where the “plus one” prevents division by 
zero. 

As the pollution rule requires that we know both the sugar consumed and sugar 
metabolised during the last move of an agent to that location it is simpler to 
incorporate the Pollution Formation rule into the movement rule. The alternative 
is to track the sugar consumed during each move which would require another 
attribute dehned in the Agent schema. 
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— Ad OVCTTlCTltpQiiii^lj^Qyi - 

ASugarScape 
step’ = step 

maxSugar' = maxSugar 
sex' = sex 

population' = population 
vision' = vision 
age' = age 
maxAge’ = maxAge 
agentCulture' = agentCulture 
loanBook' = loanBook 
children' = children 
agent Immunity' = agent Immunity 
diseases' = diseases 
metabolism' = metabolism 
initial Sugar' = initial Sugar 
Va : AGENT] I : POSITION • 

a G population' ^ distance{position'{a),position{a)) < vision{a) 

{distance{position{a),l) < vision{a) A (/ ^ run position')) 

^ [sugar{l)/{l + pollution{l))] 

< [sugar{position'{a))/{l + pollution{position'{a)))] (1) 

A {distance{l,position{a)) < distance{position'{a),position{a))) 

sugar(1 )/(I + pollution{l)) 

< sugar {position' {a)) / {1 + pollution{position' {a))) 

sugar' = sugar © {V/ ; POSITION \ I G ranposition' • / i-G- 0} 
agentSugar' = {Va ; AGENT \ a G population • 
a I—)■ agentSugar{a) + sugar {position'{a))} 

pollution' = pollution © {V/ : POSITION] x : AGENT \ position'{x) = I • 

I ^ {PRODUGTION * sugar{l) + GONSUMPTION * metabolism{x))}{2) 


(1) We use our new formula to calculate the desirability of a location; 

(2) The new pollution value for any location that an agent is present at is 
calculated as per rule dehnition. 

The rule as stated in the appendix is the generalised rule dehned for an arbitrary 
number of pollutants and resources. We have specihed the simpler version as it 
is easier to grasp. The more complex version has not been used in any of the 
Sugarscape simulations. We state the generalised rule below for completeness but 
do not present a formal specihcation of it. 
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Pollution Formation Pn,x- For n resources and m pollutants, when n- 
dimensional resource vector r is gathered from the Sugarscape the m- 
dimensional pollution production vector p is produced according to p = Hr, 
where 11 is an m x n matrix; when n-dimensional (metabolism) vector m is 
consumed then m-dimensional consumption pollution vector c is produced 
according to c = xm, where x is an m x n matrix. 


3.10. Pollution Diffusion D^- 
Pollution Diffusion Da- 

• Each a time periods and at each site, compute the pollution flux the 
average pollution level over all its von Neumann neighbouring sites; 

• Each site’s flux becomes its new pollution level. 

This rule determines how pollution diffuses over grid. Pollution diffusion is 
calculated every a turns and is computed as the average pollution level of all the 
locations von Neumann neighbours. We use the constant POLLUTIONRATE 
in place of alpha. 

The von Neumann neighbours of a location are those immediately above, below, 
left and right of the current locations (aka north, south, east and west). We dehne 
the four cardinal directions taking into account to fact that the grid wraps around 
at its edges (i.e. it is a torus). 


north : POSITION ^ POSITION 
south : POSITION ^ POSITION 
east : POSITION ^ POSITION 
west: POSITION ^ POSITION 


y x,y • 

west{{x, y)) = {{x — 1) mod M, y) 
east{{x, y)) = {{x + 1) mod M, y) 
south{{x, y)) = (x, {y — 1) mod M) 
north{{x, y)) = (x, {y + 1) mod M) 


We use this to dehne a function that returns true if two agents are von Neumann 
neighbours. It takes as parameters the two agents and a function that maps each 
agent onto their location in the simulation. 
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vonNeumanNeighbour : {AGENT x AGENT x {AGENT ^ POSITION)) boolean 

\/a,b: AGENT] position : AGENT ^ POSITION • 
vonNeumanNeighbour{a^ b, position) 
position{a) = north{position{b)) 

V position{a) = south{position{b)) 

V position{a) = east{position{b)) 

V position{a) = west{position{b)) 


_ PollutionDif fusion _ 

ALattice 

EStep 

maxSugar' = maxSugar 
sugar' = sugar 

{step mod POLLUTION RATE 7 ^ 0) =» pollution' = pollution 
{step mod POLLUTION RATE = 0) pollution' = 

{V/ : POSITION • 1 1 —)■ {pollution{north{l)) + pollution{south{l)) 
+pollution{east{l)) + pollution{west{l))) div4} 


The Pollution Diffusion Rule can be simplified slightly, by removing the re¬ 
dundant mention of flux. 

Pollution Diffusion Da- After every a time periods and at each location, 
the average pollution level over all a site’s von Neumann neighbouring 
locations becomes its new pollution level. 


3.11. Replacement - R[a^b]- 

Replacement - R[a,b] • When an agent dies it is replaced by an agent of age 
0 having random genetic attributes, random position on the Sugarscape, 
random initial endowment, and a maximum age selected from the range 
[a,b]. 

The two constants a and b we have defined already as LOWERAGELIMIT and 
UPPERAGELIMIT and we assume that the range is inclusive. It is not stated 
whether the new agents immediately consume the resources at the location they are 
placed in. We assume they do not, but accept that the alternative interpretation is 
equally valid. Although not part of the rule definition in the appendix it is stated 
elsewhere in the book that new agents will have initial resource levels set between 
5 and 25. We have defined ST ART SUG ARM IN and ST ART SUG ARM AX for 
this purpose. 

Although the simulation can be run without employing the replacement rule (in 
an effort, for example, to determine the total carrying load - maximum tolerable 
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population of agents - of a simulation space) there is no stated separate death rule. 
We will hrst add a schema that dehnes “death” explicitly to ensure consistency. 

Death : When an agent reaches its maximum allowed age or runs out of 
resources it is removed from the simulation and all its associated loans 
(either as borrower or lender) are considered void. 


_ Death _ 

A Agents 

population' = population\ 

{a : AGENT \ age{a) = maxAge{a) V agentSugar{a) = 0} (1) 

loanBook' = population' <\ loanBook > (2) 

{x : AGENT X (M X N) I first{x) G population'} 

M a-.AGENT* (3) 

a G population' 

{sex{a) = sex'{a) A vision{a) = vision'{a) 

A maxAge{a) = maxAge'{a) A agentGulture{a) = agentGulture'{a) 

A position{a) = position'{a) A age{a) = age'{a) 

A agentSugar{a) = agentSugar'{a) 

A metabolism'{a) = metabolism{a) 

A diseases' {a) = diseases{a) 

A agent Immunity'{a) = agent Immunity {a) 

A children'{a) = children{a)) 

A initial Sugar' (a) = initialSugar{a) 


(1) We remove from the population all agents who have reached their maximum 
age or who have no sugar reserves; 

(2) We remove all loans owed by or owing to these dying agents; 

(3) Any agent not being removed still has all attributes completely unchanged. 

The replacement rule follows readily from this rule, the only addition being 
the generation of new agents to replace the agents being removed. In effect we 
have broken the replacement rule into two parts Death followed by Replacement] 
although the Death rule may be used in isolation the Replacement rule must 
always be preceded by the application of the Death rule. 
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_ Replacement _ 

AAgents 

# population' = INITIALPOPULATIONSIZE (1) 

loanBook' = loanBook 
Ma:AGENT • 

a G population ^ (2) 

(a G population' 

A sex{a) = sex'{a) A vision{a) = vision'{a) 

A maxAge{a) = maxAge'{a) A agentCulture{a) = agentCulture'{a) 

A position{a) = position'{a) A age{a) = age'{a) 

A agentSugar'{a) = agentSugar{a) 

A metabolism'(a) = metabolism{a) 

A diseases' {a) = diseases{a) 

A agentlmmunity'{a) = agent Immunity {a) 

A children'{a) = children{a)) 

A initial Sugar'(a) = initial Sugar {a) 
a G population' \ population ^ 

{age'{a) = 0 (3) 

A STARTSUGARMIN < agent Sugar'{a) < STARTSUGARMAX 
A initial Sugar'(a) = agentSugar' {a) 

A diseases' {a) = 0 A children' {a) = 0) 

(1) The new population has the correct number of members; 

(2) The existing agents remain unchanged and part of the new population; 

(3) All new agents have new values initialised within the allowed limits (those 
not stated explicitly are random values within the ranges set by the speci- 
hcation invariants. 

We do not state the positions of any new agents because they are chosen randomly. 
Our schema invariants ensure that they are on the grid in a location not occupied 
by any other agent. 

We need to add some extra information to this rule dehnition to ensure that: 

(1) Newly created agents have no diseases, children or loans; 

(2) Their initial endowment of resources is within a set range. 

3.12. Agent Mating S. 

Agent Mating S: 

• Select a neighboring agent at random; 

• If the neighboring agent is of the opposite sex and if both agents are 
fertile and at least one of the agents has an empty neighboring site 
then a newborn is produced by crossing over the parents’ genetic and 
cultural characteristics; 






THE SPECIFICATION OF SUGARSCAPE 


25 


• Repeat for all neighbors. 

This rule determines how mating takes place amongst agents to produce off¬ 
spring. An agent is fertile if its age is within preset boundaries. This is represented 
by the simple isFertile function below. 


isFertile : (N x SEX) boolean 

V age : N • 
isFertile{age, male) 

MALE FERTILITY START < age < MALEFERTILITYEND 
isFertile{age, female) 

FEMALEFERTILITYSTART < age < FEMALEFERTILITYEND 

We dehne two functions that take in an agent and a mapping from parents to 
offspring and returns the father or mother of the agent. 

father ; AGENT x {{AGENT x AGENT) AGENT) AGENT 
mother : AGENT x {{AGENT x AGENT) AGENT) AGENT 

yx,mj : AGENT] Offspring : {AGENT x AGENT) AGENT • 
father{x, Of fspring) = m ^ Offspring{{m, f))=x 
mother{x, Of f spring) = f ^ Offspring{{m, f)) = x 

The issues encountered with the mating rule are similar to those with movement. 

If two sets of parent try to produce offspring in the same vacant location only one 
can succeed. As there is no preferred conflict resolution rule we cannot state any 
preference for which agents succeed in producing children and which do not. All 
we can state is that the maximum number of offspring will be produced given the 
space constraints but we cannot always be sure which offspring make it into this 
set. Neighbours in this rule refers to von Neumann Neighbours only. 

Mating although proceeding concurrently throughout the population is an ex¬ 
clusive event. That is, if agent A is mating with agent B then A cannot be mating 
with any other agent at the same time: you can only ate with one partner at a 
time. The rule itself specihes that each agent will mate with all available partners 
so the execution of the rule can involve a sequence of mating events for specihc 
agents. 

Although it is not stated in the rule dehnition the accompanying book mentions 
that each parent should gift half of its sugar to its offspring and will only mate if 
it has a sugar level equal to or greater than its initial sugar level (that is its sugar 
level on creation). This signihcantly complicates the rule and dramatically changes 
its dehnition and characteristics. However we will assume that this information 
was inadvertently omitted from the rule dehnition as the rule makes more sense if 
we include these extra factors. 
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Since each individual agent can involve itself in a sequence of up to four mating 
events during rule execution we require a specification that retains global concur¬ 
rency while still imposing a sequential ordering based on these constraints. We 
do this by collecting all possible potential mating partners into a set and then di¬ 
viding this set into a sequence of maximally sized sets where each subset contains 
only mating events that can occur concurrently. These sets are produced using a 
conflict resolution rule that ensures that only pairing that can occur simultane¬ 
ously appear within each such subset. The rule then proceeds by executing mating 
events within each subset concurrently while the sets are evaluated in sequence. 


_ Agent M ating _ 

ELattice 

AAgents 

loanBook' = loanBook 

3 potential Mating Pairs : F{AGENT x AGENT) \ 

potential Mating Pairs = {(a : AGENT,h : AGENT) \ sex{a) 7 ^ sex{b) 

A isEertile{age{a), sex{a)) A isEertile{age{head),sex{head)) 

A adjacent{position{a ), position{head ))} 

{population', position', vision', agentSugar', agentGulture', metabolism' 

, children', diseases', agentimmunity', age', sex', initialSugar') = 
concurrentM ating {getG on fictEreePair s{potentialM ating Pair s), population, position, 
agentSugar, agentGulture, metabolism, children, 
diseases, agentimmunity, age, max Age, sex, initialSugar) 


(1) Generate the set of all possible mating pairs; 

(2) Recursively proceed with concurrent mating within the conflict free subsets. 
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getConfictFreePairs : ¥{AGENT x AGENT) 

—)■ 

seq(P(AGEiVT x AGENT)) 

MAllPairs : ^{AGENT x AGENT); a : AGENT; • 

3 conflictFreeSet : F{AGENT x AGENT) \ 

conflictFreeSet C AUPairs (1) 

A a G ran conflictFreeSet a ^ domconflictFreeSet 

A a G domconflictFreeSet a ^ ran conflictFreeSet 

Mother Set : F{AGENT x AGENT) \ (2) 

other Set C AUPairs A a G doni ot/ierS'et a ^ ran AUPairs 

A a G ran other Set a ^ dom AUPairs • ^ other Set < ^ conflistFreeSet 

getGonfictFreePairs{0) = 0 

getG on f ictFreePair s{AUP airs) = (3) 

(conflictFreeSet) ^ getGonfictFreePairs{AllPairs \ conflictFreeSet) 


(1) Generate a collision free (conflict resolved) set where each agent can only 
once within the set; 

(2) Ensure this set is as large as possible; 

(3) Recurse through the remaining pairs dividing them into more conflict free 
sets. 
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concurrent Mating : sec[¥{AGENT x AGENT) 
AGENT 

X AG ENT ^ POSITION 
xAGENT 
xAGENT^N 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT ^ P AGENT 
X AG ENT -E-Pseq-B/T 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT 
xAGENT ^ SEX 
xAGENT^n 

¥ AGENT 

X AG ENT ^ POSITION 
X AGENT ^Ni 
xAGENT^N 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT ^ P AGENT 
X AG ENT -E-Pseq-BIT 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT 
xAGENT ^ SEX 
xAGENT^N 


Wtail : seqF{AGENT x AGENT)- 
head : F{AGENT x AGENT)-, 
population : P AGENT ; 
position : AGENT ^ POSITION-, 
vision -. AGENT -H- Ni; 
agentSugar : AGENT -H- N; 
agentGulture : AGENT -E- seq BIT-, 
metabolism : AGENT -H- M; 
children : AGENT ^F AGENT-, 
diseases : AGENT -H-Pseq-B/T; 
agentimmunity : AGENT -E-seq-BIT; 
age : AGENT ^ N; 
maxAge : AGENT -H- Ni; 
sex -. AGENT ^ SEX-, 
initialSugar : AGENT -H- N; 
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3newpopulation : F AGENT] 
newposition : AGENT ^ POSITION] 
newvision : AGENT Ni; 
newagentSugar : AGENT -H- N; 
newagentGulture : AGENT -H- seqBIT] 
newmetabolism : AGENT -H- M; 
newchildren : AGENT F AGENT] 
newdiseases : AGENT -H-Pseq-BIT; 
newagent Immunity : AGENT -^seqBIT] 
new age ; AGENT -H- N; 
newmaxAge : AGENT Ni; 
newsex : AGENT -H- SEX] 
newinitialSugar : AGENT -H N; | 

{newpopulation, newposition, newvision, newagentSugar, newagentGulture, 

newmetabolism, newchildren, newdiseases, new agent Immunity, newage, newmaxAge, 
newsex, newinitialSugar) = 
applyMating{asSeq{head), population, position, vision, 
agentSugar, agentGulture, metabolism, children, 
diseases, agent Immunity, age, maxAge, sex, initialSugar) • 
concurrentMating{{), population, position, vision, 
agentSugar, agentGulture, metabolism, children, 
diseases, agent Immunity, age, maxAge, sex, initialSugar) = 

{population, position, vision, agentSugar, agentGulture, metabolism, 
children, diseases, agentimmunity, age, maxAge, sex, initialSugar) 
concurrentMating{{head) ^ tail, population, position, vision, 
agentSugar, agentGulture, metabolism, children, 
diseases, agentimmunity, age, maxAge, sex, initialSugar) = 
concurrentMating{tail, newpopulation, newposition, newvision, 

newagentSugar, newagentGulture, newmetabolism, newchildren, 

newdiseases, new agent Immunity, newage, newmaxAge, newsex, newinitialSugar) 
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apply Mating : seq{AGENT x AGENT) 
AGENT 

X AG ENT ^ POSITION 
xAGENT 
xAGENT^N 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT ^ P AGENT 
X AGENT -^^FseqBIT 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT 
xAGENT ^ SEX 
xAGENT^n 

¥ AGENT 

X AG ENT ^ POSITION 
X AGENT ^Ni 
xAGENT^N 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT P AGENT 
X AG ENT -E-Pseq-BIT 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT 
xAGENT ^ SEX 
xAGENT^N 


\f population : F AGENT] 

position : AGENT ^ POSITION; 

sex : AGENT ^ SEX; 

vision : AGENT -H- Ni; 

age : AGENT N; 

initialSugar ; AGENT -H- M; 

maxAge : AGENT -H- Ni; 

metabolism : AGENT -H- M; 

agentSugar ; AGENT -E- N; 

agentGulture ; AGENT seqBIT; 

children : AGENT ^F AGENT; 

agentimmunity : AGENT seqBIT; 

diseases : AGENT -^FseqBIT; 

head : AGENT x AGENT; 

tail : seq{AGENT x AGENT); • 
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3 offspring, a, b : AGENT ; 
newsex : AGENT -H SEX] 
newvision : AGENT Ni; 

newmetabolism, newagentSugar, newinitialSugar : AGENT M; 

newmaxAge : AGENT Ni; 

newagentGulture ; AGENT -H- seqi?/T; 

newchildren : AGENT -H- ¥ AGENT] 

new a gent Immunity : AGENT -H-seq-B/T; 

inkeritedimmunity : seq BIT] 

inkeritedGulture : seqBIT] 

I offspring ^ population 
a = firstihead) A b = secondifiead) 

newckildren : ckildren U {of f spring i—)■ 0 a i—)■ ckildren{a) U [of f spring}, 
b i-A ckildrenfb) U [of f spring}} 
newsex = sex U {of f spring i— male} 

V newsex = sex U {offspring i—)■ female} 
newvision = vision U {of f spring i-A vision{a)} 

V newvision = vision U {of f spring i—)■ vision{b)} 
newmaxAge = max Age U {of f spring i-A maxAge{a)} 

V newmaxAge = max Age U {of f spring i-A maxAgefb)} 
newmetabolism = metabolism U {of f spring i—)■ metabolism{a)} 

V newmetabolism = metabolism U {of f spring i-A metabolism{b)} 
newinitialSugar = initialSugar® 

{of f spring i-A initial Sugar {a) / 2 + initialSugar{b)/2, a i—)■ initial Sugar {a)/2, b i-a initialSuga: 
newagentSugar = agentSugar U {of f spring i—)■ initialSugar} 

A Vn : 1 .. IMMUNITYLENGTH • 

{inker itedlmmunity{n) = agentlmmunity{a){n) 

V inkeritedlmmunity{n) = agentlmmunity{b){n)) 

new a gent Immunity : agentimmunity U {of f spring i-A inkeritedimmunity} 

A Vn : 1 .. GULTUREGOUNT • 

{inker itedGulture{n) = agentGulture{a){n) 

V inkeritedGulture{n) = agentGulture{b){n)) 
newagentGulture : agentGulture U {of f spring i—>• inker itedGulture} 
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applyMating{{), population, position, vision, agentSugar, agentCulture, 

metabolism, children, diseases, agentimmunity, age, maxAge, sex, initialSugar) = 
{population, position, vision, agentSugar, agentCulture, metabolism, 
children, diseases, agentimmunity, age, maxAge, sex, initialSugar) 

apply Mating {{head) ^ tail, population, position, vision, agentSugar, agentCulture, 

metabolism, children, diseases, agentimmunity, age, maxAge, sex, initialSugar) = 
if((3/oc : POSITION \ {adjacent{loc,position{ag))) 

V adjacent{loc,position{head)) A loc^ domposition) 

A {agentSugar {head) > initial Sugar {head)) A {agentSugar {ag) > initialSugar {^ 
apply Mating {tail, population U {of f spring}, position U {of f spring i—)■ loc}, 
newvision, newagentSugar, newagentCulture, new metabolism, 
newchildren, diseases U {of f spring i—)■ 0}, new agent Immunity, 
age U {of f spring i—)■ 0}, newmaxAge, newsex, initialSugar) 

else 

apply Mating {tail, population, position, vision, agentSugar, agentCulture, 

metabolism, children, diseases, agentimmunity, age, maxAge, sex, initialSugar) 

3.13. Agent Inheritance /. 

Agent Inheritance I: When an agent dies its wealth is equally distributed 
among all its living children. 

The rule definition is deceptively simple but some assumptions must be made in 
order to give it a precise definition. These assumptions are required because of the 
discrete nature of the simulation. Only living children can inherit from a parent. 

If a child is alive but scheduled to die at the same time as their parent then 
(because all agents who are due to die will die simultaneously) this child should 
not inherit from their parent. If we were to allow them to inherit we would either 
have to impose an ordering on the allocation of inheritance making the rule more 
complex or accept than the ordering will sometimes result in part of an inheritance 
disappearing. This extra complexity brings no real benefit to the simulation so we 
discount it. 

The second assumption is that we allow for rounding errors. Resources {sugar) 
come in discrete amounts so division between children requires integer division. 

This is also true of division of the loans amongst an agents children. We just 
accept any rounding errors as part of the discrete nature of the simulation. 

Finally we note that inheritance is separate from the actual death or replacement 
rule, it reallocates the resources of agents due to die but it does not remove those 
agents from the simulation. We leave that to the actual Replacement or Death 
rule and assume that one of these rules is applied after the inheritance rule. This 
simplifies the Inheritance schema. 

To enable inheritance to handle the loan book (when an agent dies its loans are 
passed on to its children) we introduce some helper functions. The asSeq function 
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turns a set of items into a sequence of items. It does not specify the ordering in 
the sequence. 

asSeq ; PX —)■ seqX 

V x : P X; y ■. seq X • 

asSeq{x) = y ^ (ran?/ = xA4^y = i^x) 


The second function disperseLoans takes in the loan book, a sequence contain¬ 
ing all the dying agents and the children of the agents and produces an updated 
loan book with the loans of the dying agents now dispersed amongst their children. 
To do this it employs a third function oneAgentLoans that takes in a single agent 
(who is marked for removal) the loans (in a sequence) held by that agent and the 
set containing its children. It outputs a new set of loans generated by dispersing 
all this agents loans amongst its children. In both cases we use sequences for the 
parameter we are recursing over as it makes the recursion easier to specify. 


disperseLoans : (F{AGENT x {AGENT x (M x N))) x seq AG ENT 
X{AGENT AGENT)) 

^F{AGENT X {AGENT x (N x M))) 

MLoans : F{AGENT x {AGENT x (N x N))); a : AGENT] 
tail : seq AGE NT; Ghildren : AGENT AGENT • 
disperseLoans{Loans, {), Ghildren) = Loans 

disperseLoans{Loans, (a) ^ tail, Ghildren) = 

disperseLoans{{{a} < Loans) U oneAgentLoans{a, asS'eg(ran({a} < Loans)), 
Ghildren{a)),tail, Ghildren) 


oneAgentLoans : AGENT x seq{AGENT x (M x N)) x F AGENT 
^F{AGENT X {AGENT x {N x N))) 

y a,borrower, inheritor ; AGENT; Ghildren : F AGENT; amt, dur, new Amt ; M; 

tail : seq{AGENT x (M x N)) • 
oneAgentLoans{a, {), Ghildren) = 0 

oneAgentLoans{a, {{borrower, {amt, dur))) ^ tail, Ghildren) = 

{x : AGENT I X G Ghildren • {x, {borrower, {amt div ^ Ghildren, dur)))} 

U oneAgentLoans{a, tail, Ghildren) 

The getMother and getFather functions simply take in an agent and the children 
set and hnds the mother (father) of the agent from this set. 
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getMother : AGENT x {AGENT AGENT) x AGENT ^ SEX ^ AGENT 
getEather : AGENT x {AGENT AGENT) x AGENT ^ SEX ^ AGENT 

y child, parent : AGENT] children : AGENT ^F AGENT • 

getM other {child, children, sex) = parent child G children{parent) A sex{parent) = fern 
getEather{child, children, sex) = parent child G children{parent) A sex{parent) = male 

An agent can inherit from at most two different agents, one male and one female. 

We use this to facilitate the specihcation by treating each sex separately. 


_ Inheritance _ 

AAgents 

population' = population A sex' = sex 

position' = position A vision' = vision 

age' = age A maxAge' = maxAge 

agentGulture' = agentGulture A children' = children 

metabolism' = metabolism A diseases' = diseases 

agent Immunity' = agent Immunity 

initialSugar' = initialSugar 

3 dying ; F AGENT] inheritEromEemale,inheritEromMale : AGENT -G N; 
yX : AGENT] 3p : AGENT \ x G population \ dying • (1) 

dominheritEromEemale = dominheritEromMale = population \ dying 
dying = {x : AGENT \ x G population A 

A {age{x) = maxAge{x) V agentSugar{x) = 0)} 
getMother{x, children, sex) ^ dying ^ (la) 

inheritEromEemale{x) = 0 
p = getMother{x, children, sex) A p G dying ^ 
inheritEromEemale{x) = 

agentSugar{p) div {population fl children{p) \ dying)) 

getEather{x, children, sex) ^ dying ^ {lb) 

inheritEromM ale{x) = 0 
p = getEather{x, children, sex) A p E dying 
inheritEromM ale{x) = 

agentSugar{p) div ^{population fl children{p) \ dying)) 

X G dying (3) 

^ agentSugar'{x) = 0 

X ^ dying (4) 

^ agentSugar'{x) = agentSugar{x) 

+inheritEromM ale{x) + inheritEromE emale{x) 
loanBook' = disperseLoans{loanBook, asSeq{dying), children) (5) 
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(1) First we construct the set of dying agents. Then using this set of dying 
agents we can construct two functions, one mapping amounts inherited 
from a female parent and one mapping amounts inherited from a male 
parent. These sets are then used to update the sugar of each agent; 

a) The function giving the amount each inheriting agent gets from its fe¬ 
male parent is constructed by hnding all healthy agents who have a dying 
mother and determining their share of their dying mother’s resources; 

b) The function listing amounts each agent gets from a male parent is 
constructed in an almost identical manner. 

(2) If an agent is dying its sugar level is set to zero (because it is being real¬ 
located to its children); 

(3) Otherwise the agents sugar level is its old level plus whatever it inherits 
from both dying parents; 

(4) Finally we update the loanBook using our disperseLoans function. 

3.14. Agent Culture K. 

Agent cultural transmission: 

• Select a neighboring agent at random; 

• Select a tag randomly; 

• If the neighbor agrees with the agent at that tag position, no change 
is made; if they disagree, the neighbor’s tag is flipped to agree with 
the agent’s tag; 

• Repeat for all neighbors. 

Group membership: Agents are dehned to be members of the Blue group 
when Os outnumber Is on their tag strings, and members of the Red group 
in the opposite case. 

Agent Culture K: Combination of the “agent cultural transmission” and 
“agent group membership” rules given immediately above. 

Group membership is dehned with the assumption that there are always an odd 
number of tags, tribe returns the affiliation of an agent based on the number of 
bits of each type in its culture sequence. If the majority of bits in a sequence are 
0 then it belongs to the blue tribe, otherwise it belongs to the red tribe. This is 
used by the culture rule. 

tribe : seq BIT —)• affiliation 
V aSeq : seq BIT • 

tribe{aSeq) = blue ^{aSeq > {0}) > ^{aSeq > {1}) 
tribe{aSeq) = red ^{aSeq > {0}) < ^{aSeq > {1}) 

flipTags is a recursive function that takes in a culture tag sequence belonging 
to an agent, a sequence of neighbouring agents and the mapping containing all 
agent’s culture tag sequences. It returns a new tag sequence generated by each 
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neighbouring agent flipping one bit chosen at random of the original agent’s tag 
sequence. It is aided in this by the function flipBit that takes in two bit sequences 
and returns a new sequence equal to the hrst bit sequence with one bit changed 
at random to match the other sequence at that position. 


flipBit : seq BIT x seq BIT seqi?/T 

y original, other, new : seq BIT • 
flipBit{original, other) = new 

^ original = other = ^ new A 
3 i : 0 .. original •y j : 0 .. original • 

^ 3 ^ new{i) = original{j)) A new{i) = other{i) 


flipTags : seq BIT x seq AGENT x {AGENT -A- seq BIT) fA seq i?JT 

yaSeq : seq BIT; ag : AGENT; tail : seq AGE NT; 

cultural Re sources : AGENT ^seqBIT • 
flipTags{aSeq, {), culturalResources) = aSeq 
flipTags{aSeq, {ag) ^ tail, culturalResources) = 

flipTags{flipBit{aSeq, culturalResources{ag)), tail, culturalResources) 


The sequence of neighbours is provided by the Gulture scheme which employs 
the asSeq function to convert a set of neighbours into a sequence. 
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_ Culture _ 

A Agents 

population' = population 

sex' = sex 

position' = position 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentSugar' = agentSugar 

children' = children 

loanBook' = loanBook 

diseases' = diseases 

metabolism' = metabolism 

agentimmunity' = agentimmunity 

initialSugar' = initialSugar 

Va : AGENT • a G population ^ agentCulture'{a) = 

flipTags{agentCulture{a), (1) 

asSeq{{b : AGENT \ adjacent{position{a),position{b))}), agentCulture) 

(1) For every agent a in the population we allow each other agent that counts 
a as a neighbour to flip one bit at random of a’s culture bit sequence. 

3.15. Combat C^- 

Agent Combat Ca- 

• Look out as far as vision permits in the four principle lattice directions; 

• Throw out all sites occupied by members of the agent’s own tribe; 

• Throw out all sites occupied by members of different tribes who are 
wealthier then the agent; 

• The reward of each remaining site is given by the resource level at 
the site plus, if it is occupied, the minimum of a and the occupant’s 
wealth; 

• Throw out all sites that are vulnerable to retaliation; 

• Select the nearest position having maximum reward and go there; 

• Gather the resources at the site plus the minimum of a and the occu¬ 
pants wealth if the site was occupied; 

• If the site was occupied then the former occupant is considered “killed” 

- permanently removed from play. 

reward is used by the combat rule and values a position based on its sugar 
content and the sugar reserves held by any agent at that position. The combat 
rule is really an extension of the movement rule where we are now allowed to move 
to locations occupied by other agents under certain predehned conditions. 
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reward ; POSITION x {POSITION M) x {AGENT POSITION) 
X {AGENT ^ M) X M 


yi : POSITION; sugar : POSITION N; agentSugar : AGENT N; 

positions : AGENT POSITION • 
if I G ran positions then 

rewar d{l, sugar, positions, agentSugar, GO MB AT LI MIT) = 

sugar{l) + min{{GOMBATLIMIT, agentSugar{positions {1))}) 

else 

rewar d{l, sugar, positions, agentSugar, GO MB AT LI MIT) = 
sugar {1) 

availMoves returns the set of all safe moves that an agent can make. 

availMoves : AGENT x {AGENT POSITION) x {POSITION M) x 
{AGENT N) X {AGENT seq BIT) x M 
POSITION 

y agent : AGENT; positions : AGENT h-> POSITION; vision : N; 
sugar ; POSITION >+^ N; agentSugar : AGENT hh- N; 
culture : AGENT >+^ seq BIT • 

availMoves{agent, positions, sugar, agentSugar, culture, vision) = 

{/ : POSITION; x : AGENT \ distance{l, positions {agent)) < vision (1) 
A positions{x) = I ^ {agentSugar{x) < agentSugar{agent) (2) 

A tribe{culture{x)) ^ tribe{culture{agent))) 

A {{distance{positions{x),l) < vision) (3) 

A tribe{culture{x)) ^ tribe{culture{agent))) ^ 
agentSugar{x) < agentSugar {agent) 

+reward{l, sugar, positions, agentSugar, GOMBATLIMIT)) • /} 

(1) Only locations within an agents neighbourhood are considered; 

(2) If a location is occupied it must be occupied by an agent belonging to a 
different tribe who has lower sugar levels; 

(3) We only consider a position already containing an agent from another tribe 
if there are no other agents from a different tribe within the neighbourhood 
of that location who are stronger than we will be once we have consumed 
the resources of the new location (that is agents who may retaliate against 
us for killing an agent belonging to their own tribe). 

We note that the rule as stated means we consider retaliation under all conditions 
even if we are just moving to an empty location. It is unclear from the definition 
given as to how exactly we check for retaliation. Do we base our check on agents 
visible from our current position or from the proposed position. We have assumed 
that it is based on the proposed position but it could easily be otherwise. We 
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also assume that the range used is based on the vision of the moving agent as this 
seems logical. 

The synchronous version of the combat rule assumes that all combat occurs 
instantaneously (concurrently). We note that it is simpler to specify in that we 
just state the before and after states and make no mention of orderings of combat. 


- C OTflbClt . 

ASugarScape 


step' = step 

maxSugar' = maxSugar 
pollution' = pollution 

loanBook' = population' < loanBook > {population' < {ranloanBook)) 
population' C population 

sugar' = sugar © {p : POSITION \ p G lanposition' • p i-A 0} 


( 1 ) 

( 2 ) 

( 3 ) 


'iag : AGENT] I : POSITION • 

ag G population' (4) 

{sex'{ag) = sex{ag) 

A vision'{ag) = vision{ag) 

A age'{ag) = age{ag) 

A max Age' {ag) = maxAge{ag) 

A children'{ag) = children{ag) 

A agentCulture'{ag) = agentCulture{ag) 

A agentimmunity'{ag) = agentlmmunity{ag) 

A metabolism'{ag) = metabolism{ag) 

A diseases' {ag) = diseases{ag) 

A initial Sugar' {a) = initialSugar{a) 

A agentSugar'{ag) = agentSugar{ag) (5) 

Preward{position'{ag), sugar, position, agentSugar, COMBATLIMIT) 

A position'{ag) G (6) 

availMoves{ag, position, sugar, agentSugar, agentCulture, vision{ag))) 

) 

ag G population \ population' ^ (7) 

3x : AGENT • position'{x) = position{ag) A tribe{culture{x)) ^ tribe{culture{ag)) 


{I G availMoves{ag,position, sugar, agentSugar, agentCulture, vision{ag)) (8) 

A reward{l, sugar, position, agentSugar, COMBATLIMIT) 

> reward{position'{ag), sugar, position, agentSugar, COMBATLIMIT) 

A distance{position{ag), 1) < distance{position{ag),position'{ag))) 

=P 3x : AGENT • position'""{1) = x A position{x) ^ I 
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(1) Every agent that is removed from the simulation is also removed from the 
loanBook] 

(2) No new agents are introduced; 

(3) Location sugar levels are updated; 

(4) Every agent that remains in the population has all its attributes unchanged 
apart from (possibly) position and sugar; 

(5) We update the sugar levels of each agent using the reward function; 

(6) Every agent has moved somewhere within their old neighbourhood; 

(7) Every agent that is no longer part of the population was removed by com¬ 
bat, that is, there is another agent (the agent that killed them) now situated 
in their old position; 

(8) If a location available to an agent and the reward of that location is better 
or equal to that agent’s new position and it was closer than that agents new 
position to its old position then it must be the case that some other agent 
has just moved to that location (otherwise we would have moved there); 

We have had to make some assumptions here. It is not stated what happens 
when there are no available moves, for example if all sites are subject to retaliation. 
We have assumed that a move is preferable to staying still and that the only time 
that an agent stays in the same position is when there are no available moves. 
That is, if every site, including our current one, is subject to retaliation then we 
do not move anywhere. A more complex interpretation would be to for an agent 
that cannot escape retaliation to attack another agent anyway and hope for the 
best but purely in the interests of simplicity we have agents remain where they 
are. 


3.16. Credit Ldr- 
Credit Ldr- 

• An agent is a potential lender if it is too old to have children, in which 
case the maximum amount it may lend is one-half of its current wealth; 

• An agent is a potential lender if it is of childbearing age and has wealth 
in excess of the amount necessary to have children, in which case the 
maximum amount it may lend is the excess wealth; 

• An agent is a potential borrower if it is of childbearing age and has 
insufficient wealth to have a child and has income (resources gathered, 
minus metabolism, minus other loan obligations) in the present period 
making it credit-worthy for a loan written at terms specihed by the 
lender; 

• If a potential borrower and a potential lender are neighbors then a 
loan is originated with a duration of d years at the rate of r percent, 
and the face value of the loan is transferred from the lender to the 
borrower; 
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• At the time of the loan due date, if the borrower has sufficient wealth 
to repay the loan then a transfer from the borrower to the lender is 
made; else the borrower is required to pay back half of its wealth and 
a new loan is originated for the remaining sum; 

• If the borrower on an active loan dies before the due date then the 
lender simply takes a loss; 

• If the lender on an active loan dies before the due date then the bor¬ 
rower is not required to pay back the loan, unless inheritance rule I is 
active, in which case the lender’s children now become the borrower’s 
creditors. 

totalOwed calculates the total amount owed from a given sequence of loans. We 
have assumed that interest is simple interest and not compound. 


totalOwed : seq{AGENT x (N x N)) ^ M 

Va : AGENT] amt,dur : N; tail : seq{AGENT x (N x N)) • 
totalOwed{{)) = 0 

totalOwed{{{a, {amt, dur))) ^ tail) = {amt + amt * RATE * DU RATION) 
+totalOwed{tail) 


canLend and willBorrow are simple rules. The dehnition of what determines 
credit-worthiness is missing so we have assumed it means an agent has enough 
money to pay all their outstanding loans. 


canLend : N x SEX x N —)■ boolean 


V age, sugar : M • 
canLend{age, male, sugar) 

age > MALEEERTILITYEND 

V {MALEEERTILITYSTART < age < MALEEERTILITYEND 

A sugar > OHILDAMT) 
canLend{age, female, sugar) 

age > EEMALEEERTILITYEND 

V {EEMALEEERTILITYSTART < age < EEMALEEERTILITYEND 

A sugar > OHILDAMT) 
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will Borrow : M x SEX x M x F{AGENT x (N x N)) —)■ boolean 

\/age, sugar : M; loans : F{AGENT x (N x N)) • 
will Borrow {age, male, sugar, loans) 4=^ 

{MALEEERTILITYSTART < age < MALEEERTILITYEND 
A sugar < GEIILDAMT) 

A sugar > totalOwed{asSeq{loans)) 
will Borrow {age, female, sugar, loans) 

{EE MALEEERTILITY ST ART < age < EEMALEEERTILITYEND 
A sugar < GHILDAMT) 

A sugar > totalOwed{asSeq{loans)) 

amt Avail depends on whether an agent can still have children. If they are no 
longer fertile then they can loan out half their available sugar. If that are still 
fertile then they have to retain enough sugar to have children. 

amt Avail : M x SEX x M N 

V age, sugar : M • 
amtAvail{age, male, sugar) = 

ii{age > MALEEERTILITYEND)then 
sugar div 2 

else a {is Eertile{age, male) A sugar > G HI LDAMT)then 
sugar - GHILDAMT 

else 

0 

amt Avail{age, female, sugar) = 
if {age > EEMALEEERTILITYEND)then 
sugar div 2 

elseif{isEertile{aqe, female) A sugar > GHILDAMT)then 
sugar - GHILDAMT 

else 

0 

amtReq is the amount that a lender requires. This is not defined so we can 
only use a best guess as to what it is. We assume that the amount required is 
that which gives the borrower enough sugar to have children. This is the simplest 
sensible dehnition we can think of. 

amtReq : N —)■ N 

V sugar : N • 

amtReq{sugar) = GHILDAMT — sugar 

We supply some simple helper functions that extract the borrower and lender 
from a loanBook entry, calculate the amount due from a loan, the principal and 
the due date (defined as the step when payment is due). 
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lender : AGENT x {AGENT x (N x M)) ^ AGENT 
borrower : AGENT x {AGENT x {N x N)) ^ AGENT 
amtDue : AGENT x {AGENT x {NxN))^N 
principal : AGENT x {AGENT x (M x M)) ^ N 
due : AGENT x {AGENT x (M x N)) ^ M 

V/,6 : AGENT] : N • 
lender{l, {b, {p, d))) = I 
borrower{l, {b, {p, d))) = b 

amtDue{l, {b, {p, d))) = p + p* RATE * DURATION 
principal{l, {b, {p,d))) = p 
due{l, {b, {p, d))) = d 

Finally, using these functions we can present the Pay Loans schema. 


_ Pay Loans _ 

AAgents 

EStep 

population' = population 

sex' = sex 

position' = position 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentGulture' = agentGulture 

a gent Immunity' = agentimmunity 

children' = children 

diseases' = diseases 

metabolism' = metabolism 

initialSugar' = initialSugar 

3dueLoans,newLoans : {AGENT FA {AGENT x (N x M))) • 

dueLoans = loanBook > {iaii{loanBook) > {a ; (M x N) | second{a) = step}) 

{loanBook', agentSugar') = 

payExclusiveLoans{chooseGonflictEreeSets{dueLoans), agentSugar, loanBook) 


This schema is complicated by the fact that it is possible that an agent has a 
loan due and cannot pay this loan off. In this case, according to the rule dehni- 
tion, the borrower must pay half of its sugar to the lender and renegotiate another 
loan to cover the remainder of its debt. Under this rule some issues will arise if 
the borrower has more than one due loan and cannot pay these loans off. The 
lender must pay each borrower in sequence the amount of half its sugar. This 
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cannot be performed simultaneously (for example if we owe three loans we can¬ 
not give each lender half our sugar as this would mean giving out more sugar 
than we actually have). In order to remain true to the rule dehnition we must, 
when we have more than one loan due, pay each loan in some sequence (dehned 
using a conflict resolution rule e.g. pay biggest loan hrst). The helper function 
chooseConflictFreeLoans returns a sequence of groups of loans that are conflict 
free (i.e. a borrower can only appear once in each group). 

The function payExclusiveLoans takes in this sequence of loan sets and pro¬ 
cesses each set concurrently in the same manner as the Mating rule. 

chooseConflictFreeLoans : {AGENT •(-)■ {AGENT x (N x N))) 
t-)- 

seq{AGENT ^ {AGENT x (N x N))) 

Va : AGENT- dueLoans : {AGENT ^ {AGENT x (M x M))) • 
chooseConflictFreeLoans{0) = {) 
chooseC on f lictFreeLoans{dueLoans) = 

3maxSet : {AGENT ■H- {AGENT x (N x N))) | maxSet C dueLoans 
A #({a} < {ran dueLoans)) > 0 ^ #({«} <1 {ranmaxSet)) = 1 (1) 

{maxSet) ^ chooseConflictFreeLoans{dueLoans \ maxSet) 

(1) We choose the largest convict free set possible where a set is deemed conflict 
free if all borrowers only appear in that set at most once. 


payExclusiveLoans : seq{AGENT -H- {AGENT x (M x N))) 
xAGENT^n 

x{AGENT^ {AGENT x (N x M))) 

t-)- 

{{AGENT o {AGENT x (N x M))) 

X AGENT ^N) 

Vtail : seq{AGENT ^ {AGENT x (M x M))); 
head,loanBook ; {AGENT x {AGENT x (N x N))); 
agentSugar ; AGENT -P- N • 
payExclusiveLoans{{), agentSugar, loanBook) = 

{loanBook, agentSugar) 

payExclusiveLoans {{head) ^ tail, agentSugar, loanBook) = 

3 new Agent Sugar : AGENT -P- N; newLoans : {AGENT <H- {AGENT x (M x M))) 
{newLoans,new AgentSugar) = makePayments{asSeq{head), 0, agentSugar) • 
pay ExclusiveLoans{tail, new Agent Sugar, {LoanBook \ head) U newLoans) 

makePayments is a recursive function that goes through a sequence of loans 
and makes the hnal payment on each one. It is used in the PayLoans schema 
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where it takes in a sequence of the due loans and the agents current sugar levels 
and returns a set of renegotiated loans, where payment is unable to be made, and 
the new agent sugar levels. 


makePayments : sec\{AGENT x (AGENT x (N x N)))x 
F{AGENT X {AGENT x {N x N))) x {AGENT ^ N) 

-)■ 

{^{AGENT X {AGENT x (N x M))) x {AGENT M)) 

yrenegotiatedLoans,new : F{AGENT x {AGENT x (N x M))); 
updatedSugar, agentSugar : AGENT -H N; 
loan : {AGENT x {AGENT x (M x M))); 
tail : seq{AGENT x {AGENT x (N x M))) • 

makePayments{{),renegotiatedLoans,updatedSugar) = (1) 

{renegotiatedLoans, updatedSugar) 

makePayments{{loan) tail, new, agentSugar) = (2) 

if amtDue{loan) < agentSugar{borrower{loan)) then (2a) 

makePayments{tail, new, agentSugar 

® {lender{loan) i—)■ agentSugar{lender{loan)) + amtDue{loan), 
borrower {loan) i—)■ agentSugar{borrower{loan)) — amtDue{loan)}) 
else (26) 

makePayments {tail, new U {{lender {loan), 

{borrower{loan), {amtDue{loan) — agentSugar{borrower{loan)) div 2, 
due{loan) + DURATION)))}, 
agentSugar © {lender{loan) i—)■ agentSugar{lender{loan)) 

PagentSugar{borrower{loan)) div 2, 

borr ower {loan) ha agentSugar{borrower{loan)) div 2}) 

For the hnal part of the Gredit rule we need to be able to work out the total 
owed by an agent over all loans. First we dehne two helper functions: sumLoans 
and totalOwed. 


sumLoans : seq{AGENT x {AGENT x (N x N))) 

ytail : seq{AGENT x {AGENT x (M x M))); 

top : AGENT x {AGENT x (N x N)) • 
sumLoans{{)) = 0 

sumLoans {{top) ^ tail) = sumLoans {tail) + amtDue{top) 
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totalOwed ; AGENT x {AGENT x {AGENT x (N x N))) ^ M 
total Loaned : AGENT x {AGENT x {AGENT x {N x N))) ^ 

y agent : AGENT- loans : AGENT x {AGENT x (N x M)) • 
totalOwed{agent, loans) = 

sumLoans{asSeq{loans > ({a^fent} < {ranloans)))) 
totalLoaned{agent, loans) = 

sumLoans{asSeq{{agent} <\ loans)) 
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_ M akeLoans _ 

A Agents 
EStep 

population' = population A sex' = sex 
position' = position A vision' = vision 
age' = age A maxAge' = maxAge 
initialSugar' = initialSugar 

agentCulture' = agentCulture A agentimmunity' = agent Immunity 
diseases' = diseases A children' = children A metabolism' = metabolism 


3 new Loans : F{AGENT x {AGENT x (N x N))); 

V a^f, lender, borrower : AGENT ; amt, due : N • 

loanBook' = loanBook U newLoans (1) 

ag G doTunewLoans 

agentSugar'{ag) = agentSugar{ag) — totalLoaned{ag, newLoans) (2a) 

G dom(ran newLoans) {2b) 

agentSugar'{ag) = agentSugar{ag) + totalOwed{ag, newLoans) 
ag ^ dom{newLoans) U dom(ran newLoans) (2c) 

agentSugar' {ag) = agentSugar{ag) 
willBorrow{age{ag), sex{ag), agentSugar'{ag), 

ran(Zoan5ooA:' n {a : AGENT x {AGENT x (M x M)) 

I borrower {a) = borrower {loan)})) {2d) 


-1 3ag2 : AGENT • canLend{age{ag2), sex{ag2), agentSugar'{ag2)) 
A adjacent{position{ag2),position{ag)) 

totalLoaned{ag, newLoans) < amtAvail{age{ag), sex{ag), agentSugar{ag)){3) 


totalOwed{ag, new Loans) < amtReq{agentSugar{ag)) (4) 

{lender, {borrower, {amt, due))) G newLoans (5) 

{canLend{age{lender), sex{lender), agentSugar{lender)) (5a) 

A willBorrow{age{borrower), sex{bor rower), agentSugar {borrower), 
{borrower} < {ran loanBook)) 

A amt < min{{amt Avail{age{lender), sex{lender), agentSugar {lender)), 

amt Req{agent Sugar {borrower))}) {5b) 

A due = step + DURATION (5c) 

A adjacent{position{lender),position{borrower))) {5d) 


(1) The new loan book is the old book plus the new loans; 

(2) The following properties ensure sugar is updated correctly and that the 
correct amount of borrowing has taken place: 

a) If an agent is a lender then their new sugar levels decrease by the amount 
the have lent; 

b) If an agent is a borrower then their sugar has increased by the amount 
they have borrowed; 
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c) Any agent that neither borrowed or lent has the same sugar levels as 
before; 

d) If there remain any agents who still need to borrow then it is because 
there are no agents in their neighbourhood who are still in a position to 
borrow. 

(3) The total amount loaned by any agent is no greater than the amount that 

agent had available; 

(4) The total amount borrowed is less than or equal to the amount required 

by the borrower; 

(5) Every loan in this set must have the following properties: 

a) The lender must be in a position to lend; 

b) The borrower must need to borrow; 

c) The amount is less than or equal to the minimum of (i) the amount 
required by the borrower and (ii) the maximum amount available from 
the lender; 

d) The due date of the loan is set by the DURATION constant; 

e) the borrower and lender must be neighbours. 


3.17. Agent Disease E. 

Agent immune response: 

• If the disease is a substring of the immune system then end (the agent 
is immune), else (the agent is infected) go to the following step; 

• The substring in the agent immune system having the smallest Ham¬ 
ming distance from the disease is selected and the first bit at which it 
is different from the disease string is changed to match the disease. 

Disease transmission: For each neighbor, a disease that currently afflicts 
the agent is selected at random and given to the neighbor. 

Agent disease processes E: Combination of “agent immune response” and 
“agent disease transmission” rules given immediately above 

subseq is a function for determining whether one sequence is a subsequence of 
another, hammingDist determines the number of bit differences in two sequences 
of the same size. 


subseq : seq BIT x seqBIT —)■ boolean 

\/mid, aSequence : seqBIT • 
subseqimid, aSequence) 

3prefix, suffix : seqBIT • prefix ^ mid ^ suffix = aSequence 
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hammingDist : seq BIT x seqBIT —N 

Wtail, rest : seqBIT • 
hamming Dist {{), ()) = 0 

hamming Dist {{1) ^tail, (1) rest) = hamming Dist{tail, rest) 

hamming Dist {{0) ^ tail, (0) rest) = hamming Dist{tail, rest) 

hamming Dist {{t)) ^ tail, (1) ^ rest) = 1 + hamming Dist {tail, rest) 

hamming Dist {{1) ^ tail, (0) ^ rest) = 1 + hamming Dist {tail, rest) 

apply Diseases takes in a bit sequence representing the immunity of an agent and 
a list of diseases that affect the agent and produces a new immunity bit sequence 
that is updated by the disease list. More precisely, for every disease not in the 
immunity sequence a single bit in the closest subsequence that matches the disease 
is flipped to make the sequence more closely match the disease. It uses another 
function processinfection to process each disease in the disease set. 


apply Diseases : seqBIT x seq seq BIT —>■ seqBIT 

'i I, d\ seqBIT] tail : seq seq BIT • 
applyDiseases{I, {)) = I 
applyDiseases{I, {d) ^ tail) = 

applyDiseases{processInf ection{I, d),tail) 


processinfection : seq BIT x seqBIT seqBIT 

yi,d: seqBIT • 
if subseq{d, /)then 

processInfection{I, d) = I 

else 

3a,b,c : seq BIT ; V x : seq BIT • 
a^h^c=I 

(#^ = #3; A subseq{x, I)) hamming Dist{b,d) < hamming Dist {x,d) 

: 1 .. # / • {y{'i) ^ d{'i) A Vj : N • j < i ^ d{j) = y{j)) 
processInfection{I, d) = I © {(i + # a) h-)■ b{i)}) 


ImmuneResponse is the simplest part of this rule to specify. The recursive 
function applyDiseases does all the work. 
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_ Immune Re spouse _ 

AAgents 

loanBook' = loanBook 
population' = population 
sex' = sex 
position' = position 
vision' = vision 
age' = age 
maxAge' = maxAge 
agentCulture' = agentCulture 
diseases' = diseases 
children' = children 
agentSugar' = agentSugar 
metabolism' = metabolism 
initial Sugar' = initial Sugar 

agent Immunity' = {a : AGENT \ a G population • 

a i-G- applyDiseases{agentlmmunity{a), asSeq{diseases(a)))} 

\/X : AGENT • x E population ^ agentSugar'{x) = agentSugar{x) — 

ih{d : seq-B/T | d G diseases{a) A subseq{d, agent Immunity {a)} (1) 


Although not stated in the rule definition careful reading of the accompanying 
text [Epstein and Axtell, 1996] shows that there is a penalty that is applied to 
each agent carrying diseases that it has no immunity to. The text states that for 
every disease carried by an agent that it has no immunity to, sugar metabolism 
is increased by one. So if an agent carried two diseases that it has no immu¬ 
nity to then its metabolism rate increases by two. This extra cost can equally 
be deducted by the metabolism rule or the disease rule. Purely for the sake of 
narrative it is placed in the ImmuneResponse rule where it is first referenced in 
the original Sugarscape book. This is implemented by the final two lines (1) of the 
ImmuneResponse schema. 

The transmission of diseases is the more complex part of this rule. We will 
again use a recursive helper function newDiseases to construct a set of diseases 
that an agent can catch from its neighbours. It takes the set of neighbours and 
their current diseases as input and constructs a set of diseases where one disease 
is chosen from each neighbour. 
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new Diseases : seq AGE NT x {AGENT -H F{seq BIT)) —)■ Fseq BIT 

Va : AGENT] tail : seq AGENT] diseases : AGENT -H-P(seq BIT) • 
newDiseases{{), diseases) = 0 
new Diseases {{a) ^ tail, diseases) = 
if diseases{a) = 0then 

newDiseasesitail, diseases) 

else 

3d : seqBIT \ d G diseases{a) • {d} U newDiseasesitail, diseases)) 


_ Transmission _ 

A Agents 

loanBook' = loanBook 
population' = population 
sex' = sex 
position' = position 
vision' = vision 
age' = age 
maxAge' = maxAge 
agentGulture' = agentGulture 
agentimmunity' = agentimmunity 
children' = children 
agentSugar' = agentSugar 
metabolism' = metabolism 
initialSugar' = initialSugar 
V a : AGENT • a G population 

diseases' (a) = diseasesia) U (1) 

newDiseasesiasSeqivisibleAgentsia,position, 1)), diseases) 


(1) visibleAgents returns the set of neighbours of an agent and this set is then 
passed to the newDisease function which returns a set of diseases, one 
chosen from each agent in the neighbour set. 


3.18. Rule Application Sequence. Each tick of simulation time consists of the 
application of a set sequence of rules. Not all rules can be used together so we 
identify the allowable sequences of rules. We note that it is not stated in the book 
what order the rules are to be applied. In the absence of this information we will 
pick one ordering and restrict ourselves to this ordering. 

The display the different allowable combinations of rules in any given simulation 
we use the following terminology. 
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{Rule}: The indicates that Rule is optional. We can choose to include it or 
not in a simulation; 

Rule A I RuleB: This indicates that there is a choice of which rule to apply 
- either one or the other but not both. 

This rule ordering is for simulations using only a single resource and so omits the 
Trade rule. 


Tick 

[9 Growback I9 SeasonalGrowback] 

[9 Movementbasic \ (9 MovementpoUuUon 9 PollutionDiffusion) I9 Gombat] 

{9 Inheritance}}^ Death{[^ Replacement I9 AgentMating]}} 

{9 Gulture}}^ Pay Loans 9 MakeLoans} 

{9 Transmission 9 ImmuneResponse} 

4. Asynchronous Sugarscape Specification 

AU is the sequential application of rules to agents during a simulation step. 
If, for example, all agents move during a single step then a sequential ordering is 
imposed on all of the agents and they will move one at a time (that is, sequentially) 
based on that ordering. This is in contrast to SU where all agents will attempt 
to move simultaneously (concurrently). AU is easier to implement that SU as it 
maps directly onto the current standard sequential programming practice. AU 
requires no collision detection and resolution (as for example when two agents try 
to simultaneously move to the same location) because concurrency is excluded - 
only one agent can move at any one time. It is well know that the AU and SU 
approaches can deliver different simulation results. 

Although AU and SU are both commonly used in CA based simulations, where 
agent interactions are simple in nature, AU is prevelent in ABM. This is due to 
the lack of any good SU algorithms that can handle the complex interactions such 
as Movement, Gombat or Trade that can appear in ABM based simulations. 

We have provided a specihcation of Sugarscape that assumes SU. For the sake of 
completeness and to allow us to make comparisons between synchronous and asyn¬ 
chronous updating in Sugarscape we will now present an AU based specihcation 
of the rules of Sugarscape. 

4.1. Variants of Asynchronous Updating. There are a number of varieties of 
AU. These variations differ in how they sequentially order agents for updating. 
The best known variations are [?]: 

Fixed Direction Line-By-Line: The locations in the lattice representing 
the simulation space are updated in the order they appear in the lattice 
(usually left to right, top-down); 
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Fixed Random Sweep: The order that is used is determined randomly 
at the start of the simulation and this order is used for every step in the 
simulation; 

Random New Sweep: The order that the agents are updated in is deter¬ 
mined randomly at the start of each step (each step uses a different order); 

Uniform Choice: Each agent has an equal probability of being chosen. 

If there are n agents, then n agents are chosen randomly during a step. 

During any single step an agent may not be picked at all or may be picked 
more than once (in contrast Random New Sweep guarantees every agent is 
picked exactly once per step); 

Exponential Waiting Time: This is a Time Driven method, all the oth¬ 
ers are step driven. Every agent has its own clock which rings when the 
agent is to be updated. The waiting times for the clock are exponentially 
distributed (with mean 1). The probability that an event occurs at time t 
follows where f is a real number, f > 0. This is most similar to Uniform 
Choice. 

We will provide a specihcation for each variation in turn. 

Fixed Direction Line-by-Line takes in a set of agents and their positions on the 
lattice. It produces a sequence of agents where every agent appears once and only 
once in the sequence and the order of the sequence is determined by the agents 
position on the lattice. 

lineByLine : AGENT ^ POSITION 
^seqAGENT 

\f thePositions : AGENT POSITION] theSequence : seq AGENT • 
lineByLineitheSet) = theSequence 

T an the Sequence = dovatheP ositions A theSequence = theP ositions (1) 

(n, a) G theSequence n = first(thePositions{a)) * DIM + secondithePositions{a)) (2) 

(1) Each agent in the population appears in the sequence once and only once; 

(2) If one agent appears before another in the sequence then it also appears 
before that agent on the lattice. 

Fixed Random Sweep returns a sequence of the agents in some hxed random 
ordering. This random ordering is chosen once at the start of the simulation and 
is hxed for the entire simulation run. 

RANDOMORDER : seqPOSITION 

# RANDOMORDER = # POSITION (1) 

Vn, m ; N • RANDOMORDER{n) = RANDOMORDER{m) (2) 

^ n = m 
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(1) RANDOMORDER is a globally defined sequence containing an ordering 
of positions on the lattice; 

(2) Each position on the lattice appears once and only once in this sequence. 

Any ordering that satisfies these constrains is allowable according to our specifi¬ 
cation. This introduces the randomness into the sequence. 


fixedRandom : AGENT ^ POSITION 
^seqAGENT 

\/ thePositions : AGENT ^ POSITION; theSequence : seq AGENT • 
fixedPandomithePositions) = theSequence 

I an the Sequence = domtheP ositions A ^theSequence = j^theP ositions (1) 
Vi : 0 .. # theSequence — 2; Oi, 02 : AGENT • 

(i,ai) G theSequence A (i -f l,ai) G theSequence 
(3xi,X 2 : N I (xi, ai), (x 2 , 02 ) G RAN DOMORDER 

A xi < X2 (2) 

(1) Every agent in the population appears once and only once in the resulting 
sequence; 

(2) The ordering of agents in the sequence is based on the ordering defined in 
RANDOMORDERING. 

Random New sweep is simpler to specify. We return a random ordering of agents 
after each call. We only need to ensure that every agent appears in this sequence 
exactly once. 


rndNewSweep : AGENT ^ POSITION 
^seqAGENT 

y thePositions : AGENT ^ POSITION; theSequence : seq AGENT • 
rndN ew SweepitheP ositions) = theSequence 

ran theSequence = dovatheP ositions A 4!^ theSequence = 44 thePositions{l) 

(1) Every agent in the population appears once and only once in the resulting 
sequence; 

Uniform Choice allows for an agent to be picked multiple times. The only 
constraints are that the sequence returned contains only agents in the population 
and that the size of the sequence equals the number of agents. 
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uniformChoice : AGENT ^ POSITION 
•H- seq AGENT 

'i theP ositions : AGENT ^ POSITION] n : N; theSequence : seq AGENT | 0 < n < jj^theSequ 
uniformGhoiceithePositions) = theSequence 
theSequence{n) G domthePositions 
A the Sequence = ij^theP ositions 

(1) Every agent in the sequence is an agent from the simulation population; 

(2) the size of the sequence equals the total number of individual agents in the 
population. 

Each variation of asynchronous updating can now be covered by the simple 
matter of swapping in the appropriate ordering function within the specihcations. 


4.2. Growback, Seasonal Growback and Replacement. Replacement, Growback 
and SeasonalGrowback belong to the category of rules we term independent. This 
category includes all rules where the agent involved in the update (or rule execu¬ 
tion) does not interact with any other agent - the update result is independent of 
any outside factor. It follows then that the order in which these rules are executed 
will have no bearing on their outcome. Given this we need make no changes to 
any of these rules. 


4.3. Pollntion Diffusion. PollutionDiffusion is dehned specihcally as a syn¬ 
chronous rule. There is no asynchronous alternative to this rule as imposing AU 
would redehne the rule entirely. For this reason we do not produce a AU specih- 
cation of this rule. 


4.4. Movement. The specihcation of rules under an AU regime follows a standard 
pattern. First we impose an ordering on all the agents subject to the rule and then 
we recursively apply the update to each agent in the dehned order. Each individual 
agent update can affect the global state and these changes must be passed forward 
to the next sequence of agent updates. This is in contrast to SU where all updates 
occur simultaneously. 

We always dehne the application of the rule to agents in a sequence recursively. 
While the rules themselves can be quite simple the Z notation forces us to pass to 
each update all parts of the global state that can be changed. This can result in 
large function signatures. 
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_ AsyncM ovementi^asic _ 

ASugarScape 

step' = step 
lod = loc 

maxSugar' = maxSugar 

pollution' = pollution 

sex' = sex 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentCulture' = agentCulture 

loanBook' = loanBook 

diseases' = diseases 

agent Immunity' = agent Immunity 

children' = children 

metabolism' = metabolism 

population' = population 

initial Sugar' = initial Sugar 

{sugar', agentSugar', position') = 

applyMove{rndNewSweep{position), vision, sugar, agentSugar, position) 


Movement is a typical example of this structure. The main specihcation AsyncMovementb, 
simply passes the relevant state information alongside the ordering of agents (ac¬ 
cording to whatever AU variant we are using) to the recursive function apply Move. 

This recursive function applies the move rule to each agent in turn and returns 
the hnal updated agent position, agent sugar levels and lattice sugar levels. 
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apply Move : seq AGENT 
xAGENT^N 
x{POSITION y+^N) 
xAGENT 

X AGE NT POSITION 

FA 

{{POSITION 
xAGENT y+^N 
X AGE NT POSITION) 

y head : AGENT] tail : seq AGENT] population : ¥ AGENT] 
positions : AGENT POSITION] sugar ; POSITION 
agentSugar : AGENT N; vision : AGENT -A N; • 

applyMove{{), vision, sugar, agentSugar,positions) = (1) 

{sugar, agentSugar, positions) 

applyMove{{head) ^ tail, vision, sugar, agentSugar, positions) = (2) 

3newLoc : POSITION \ newLoc G neighbourhood{position{head),vision{head)) (3) 
A y otherLoc : POSITION \ otherLoc G neighbourhood{position{head), vision{head)) 
sugar {other Loc) < sugar {newLoc) • 
applyMove{tail, vision, sugar © {newLoc i-A 0 }, 

agentSugar © {head i—)■ agentSugar{head) + sugar {new Loc)}, 
positions © {head i—)■ new Loc}) 


(1) The base case. If there are no agents left to update then we simply return 
the current state; 

(2) The recursive case. If we have agents left to process then we move the hrst 
agent in the list and apply the rule to the remaining agents; 

(3) Find the best location for the agent to move to based on sugar levels at 
each location. 


4.5. Pollution Diffusion. The movement rule for pollution is almost identical 
to the simpler basic movement rule. It only differs in that it takes pollution into 
account when selecting the best new position for an agent to move to. 
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AsyncM ovementpoiiuUon - 

ASugarScape 

step' = step 
lod = loc 

maxSugar' = maxSugar 
sex' = sex 

pollution' = pollution 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentCulture' = agentCulture 

loanBook' = loanBook 

children' = children 

agent Immunity' = agent Immunity 

diseases' = diseases 

metabolism' = metabolism 

population' = population 

{sugar', agentSugar', position') = (1) 

applyMovCpoiiutioni'f'ndNewSweep{position), vision, sugar, 
agentSugar, position, pollution) 

pollution' = pollution® {V/ : POSITION] x : AGENT \ position'{x) = / • (2) 
I ^ {PRODUCTION * sugar{l) + CONSUMPTION * metabolism{x))} 


(1) Call the recursive applyMovepoUuUon to apply movement rule to each agent 
in turn; 

(2) Update location pollution levels based on agent movement. 
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apply Movepoiiution ■ seq AG ENT 
xAGENT^N 
x{POSITION y+^N) 
xAGENT 

X AGE NT POSITION 

xPOSITION>+^N 

FA 

{{POSITION >+^N) 
xAGENT y+^N 
X AGE NT POSITION) 

V head : AGENT ; tail : seq AGENT ; population : P AGENT ; 
positions : AGENT POSITION] sugar ; POSITION N; 
agentSugar : AGENT y+-y N; 

vision : AGENT -H N; pollution : POSITION h-> N • 
applyMovepoiiution{{)^ vision^ sugar, agentSugar, positions, pollution) = 
{sugar, agentSugar, positions) 


applyMovepoiiution{{head) ^ tail, vision, sugar, agentSugar, positions, pollution) = 
3newLoc : POSITION \ newLoc G neighbourhood{position{head),vision{head)) 

A y otherLoc : POSITION \ otherLoc G neighbourhood{position{head), vision{hea 
sugar {other Loc) div (1 + pollution{other Loc)) < 
sugar {new Loc) div (1 + pollution{position' {newLoc))) • 
applyMovepoiiution{tail, vision, sugar © {newLoc i—)■ 0}, 

agentSugar © {head i—)■ agentSugar{head) + sugar {new Loc)}, 
positions © {head i-A- new Loc], pollution) 


4.6. Combat. Asynchronous Combat is undertaken with the apply AUG ombat 
function which applies the combat rule to each agent in a random order using 
the singlePight function. We note in passing that the synchronous specihcation 
seems to us to be simpler than the asynchronous one (even if the implementation 
is not). 
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_ Combat, 


async - 


ASugarScape 


step' = step 

maxSugar' = maxSugar 
pollution' = pollution 

loanBook' = population' < loanBook > {population' < {lanloanBook)) 


\/ag : AGENT] I : POSITION • 
ag G population' ^ 

{sex'{ag) = sex{ag) 

A vision'{ag) = vision{ag) 

A age'{ag) = age{ag) 

A maxAge'{ag) = maxAge{ag) 

A children'{ag) = children{ag) 

A agentCulture'{ag) = agentCulture{ag) 

A agentimmunity'{ag) = agentlmmunity{ag) 

A metabolism'{ag) = metabolism{ag) 

A diseases' {ag) = diseases{ag)) 

A initial Sugar' {ag) = initial Sugar{ag) 

{population'^ position', sugar', agentSugar') = 

applyAllCombat{rndNewSweep{position), population, position, sugar, 
agentSugar, vision, agentCulture) 
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applyAllCombat : seq AGENT 
xFAGENT 

X {AGENT POSITION) 
x {POSITION >^N) 

X {AGENT N) 

X {AGENT ^ N) 

X {AGENT ^ seq BIT) 

-)■ 

{F AGENT 

X {AGENT POSITION) 

x{POSITION>+^N) 
x{AGENT >+^N)) 

y head : AGENT] tail : seq AGENT] population : F AGENT] 

positions : AGENT POSITION] sugar ; POSITION N; agentSugar : AGENT 
vision : AGENT N; culture : AGENT seqBIT • 

apply AUG onibat{{), population, positions, sugar, agentSugar, vision, culture) = 
{population, positions, sugar, agentSugar, vision, culture) 


apply AUG ombat{{head) ^ tail, population, positions, sugar, agentSugar, vision, culture) = 
if (/lead G population)t\ien 
apply AUG ombat {tail, 

singleEight{head,population, positions, sugar, agentSugar, vision, culture)) 

else 

apply AUG ombat {tail, population, positions, sugar, agentSugar, vision, culture) 
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singleFight ; AGENT y. 

¥ AGENT 

X {AGENT POSITION) 

X {POSITION N) 

X {AGENT N) 

X {AGENT ^ N) 

X {AGENT ^seq BIT) 

-)■ 

{F AGENT 

X {AGENT POSITION) 

X {POSITION N) 
x(AGEArTH^N)) 

y agent : AGENT] population : PAG-EiVT; positions : AGENT h-> POSITION] 
sugar : POSITION NagentSugar : AGENT h-^ N; vision : AGENT -E- N; 
culture : AGENT -H- seq BIT • 

singleFight{agent, population, positions, sugar, agentSugar, vision, culture) = 
if {availMoves{agent,positions, sugar, agentSugar, culture, vision{agent)) = 0) then 
{population, positions, sugar, agentSugar) 

else 

3/oc : POSITION] available : F POSITION] 

'SotherLoc : POSITION \ 

loc, otherLoc G availMoves{agent, positions, sugar, agentSugar, culture, vision{agen 
A other Location 7 ^ location • 

reward{loc, sugar, position, agentSugar, OOMBATLIMIT) 

> reward{otherLoc, sugar, position, agentSugar, OOMBATLIMIT) 
{distance{position{agent),loc) > distance{position{otherLoc),position'{agent)) 

reward{loc, sugar, position, agentSugar, OOMBATLIMIT) 

> reward{otherLoc, sugar, position, agentSugar, OOMBATLIMIT)) • 
{population \ {positions {loc)}, {positions ^ {loc}) © {agent 1 —)■ loc}, 
sugar © {loc 1 —)■ 0}, agentSugar © {agent 1 —)■ agentSugar{agent) + 

reward{position'{agent), sugar,position, agentSugar, OOMBATLIMIT)}) 


4.7. Disease. Disease is a simple rule that follows the standard pattern for AU 
specihcation. We place all agents into a sequence, ordered according to the vari¬ 
ation of AU we are using, and apply the rule to each agent in turn updating the 
state as we go along. 




THE SPECIFICATION OF SUGARSCAPE 


63 


_ Transmission _ 

A Agents 

loanBook' = loanBook 
population' = population 
sex' = sex 
position' = position 
vision' = vision 
age' = age 
maxAge' = maxAge 
agentCulture' = agentCulture 
agentimmunity' = agentimmunity 
children' = children 
agentSugar' = agentSugar 
metabolism' = metabolism 
initialSugar' = initialSugar 
diseases' = 

applyTransmission{rndNewSweep{position), diseases, position) (1) 


(1) Call recursive applyTransmission on each agent in population in deter¬ 
mined order. 


applyTransmission : seq AGENT 
y. AG ENT P seq 5JT 

X AGE NT ^ POSITION 

FA 

AGENT ^FseqBIT 

y head ; AGENT] tail : seq AGENT] diseases : AGENT seqBIT] 
position : AGENT ^ POSITION 

3 new Infections : Pseqi?/T | newinfections (1) 

= newDiseases{asSeq{visibleAgents{head,position, 1)),diseases) • 
applyTransmission{{), diseases,position) = (2) 

diseases 

applyTransmission{{head) ^ tail, diseases, position) = (3) 

applyTransmissionitail, diseases © {head i-A- {diseases{head) U newinfections)} 
, position) 

(1) Construct a set of new infections for an agent using the previously dehned 
newDiseases] 

(2) Base case: Noting to do, return new disease mapping; 
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(3) Recursive case: Add new diseases to the first agent in the list (according 
to the rule dehnition) and then recursively apply the rule to the rest of the 

list. 

4.8. Culture. Culture is specihed in an identical manner to Disease. 

_ AsyncCulture _ 

AAgents 

population' = population 

sex' = sex 

position' = position 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentSugar' = agentSugar 

children' = children 

loanBook' = loanBook 

diseases' = diseases 

metabolism' = metabolism 

agentimmunity' = agentimmunity 

initial Sugar' = initial Sugar 

agentCulture' = applyCulture{rndNewSweep{j)osition), agentCulture^ position) 


(1) Call recursive apply Culture on each agent in population in determined 
order. 


applyCulture : seqACENT 
X AGENT ^seqBIT 
X AG ENT ^ POSITION 

AGENT ^seqBIT 

head : AGENT] tail : seq AGENT] culture : AGENT -P-seqR/T; 
position : AGENT ^ POSITION] 

3n : F AGENT \ n = neighbours{head, position{head), 1) • 
applyCulture{{), culture,position) = (1) 

culture 

applyCulture{{head) ^ tail, culture, position) = (2) 

apply Culture{tail, 

culture © {head flipTags{culture{head), asSeqin), culture)} 

, position) 

(1) Base case: return new values for culture tags; 
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(2) Recursive case: Flip the tags of the hrst agent in the list and repeat (re¬ 
cursively) for the remaining agents in the list (sequence). 


4.9. Inheritance. Inheritance also follows the same pattern as Culture and 
Disease. 


_ Asyncinheritance _ 

A Agents 

population' = population 
sex' = sex 
position' = position 
vision' = vision 
age' = age 
maxAge' = maxAge 
agentCulture' = agentCulture 
children' = children 
metabolism' = metabolism 
diseases' = diseases 
agentimmunity' = agentimmunity 
initialSugar' = initialSugar 

{loanBook', agentSugar') = (1) 

applyInheritance{rndNewSweep{position), children, loanBook, agentSugar) 


(1) Use recursive apply Inheritance function to calculate inheritance based on 
each agent in turn. 
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apply Inheritance : seq AGE NT 

X AGENT ^ P AGENT 

x{AGENT^ {AGENT x (N x M))) 

xAGENT^N 

xAGENT^N 

xAGENT^N 

-H- 

{AGENT ^ {AGENT x (N x N))) 
xAGENT^N 


yhead : AGENT] tail : seq AGENT] children : AGENT AGENT] 

loans : AGENT o {AGENT x (N x M)); 

agentSugar : AGENT -H- N; age, max Age : AGENT N 

3new Loans : AGENT -H- {AGENT x (M x N)); new AgentSugar ; AGENT -H- M 


I new AgentSugar = agentSugar 0 {{head i—)■ 0} U (1) 

{{a, amt) \ a G children{head) A amt = agentSugar{a) 
+agentSugar{head)/ # children{head)}) 

A newLoans = ({/learf} ^ loans) U oneAgentLoans{a, asS'eg(ran({a} < loans)), 

Ghildren{head)) • ( 2 ) 

applylnheritance{{), children, loans, agentSugar, age, maxAge) = (3) 

{loans, agentSugar) 

applylnheritance{{head) ^ tail, children, loans, agentSugar, age, maxAge) = 
if {age{head) = maxAge{head) V agentSugar {head) = 0) then (4) 


apply Inheritance{tail, children, newLoans, new AgentSugar, age, maxAge) 

else 

applylnheritance{tail, children, loans, agentSugar, age, maxAge) 

(1) Distribute the dying agents sugar equally amongst its children; 

(2) Distribute any loans where the dying agent is the lender equally amongst 
its children; 

(3) Base case of recursion. Nothing to do but return results; 

(4) Recursive case: If the hrst agent in the list is dying then handle that agents 
inheritance and recurse through the rest of the agents otherwise just ignore 
it and apply the rule to rest of agents. 


4.10. Mating. The AU specihcation of Mating is simpler than the SU version as 
it does not have to construct conflict free sets. It just puts all of the potential pairs 
in a sequence ordered according to the variant of AU we are using and applies the 
rule to each in turn. 
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_ AsyncAgentM ating _ 

ELattice 

^Agents 

loanBook' = loanBook 

3 potential Mating Pairs : F{AGENT x AGENT) \ 

potential Mating Pairs = {(a ; AGENT,h : AGENT) \ sex{a) ^ sex{h) 

A isEertile{age{a), sex{a)) A isEertile{age{head), sex{head)) 

A adjacent{position{a),position{head))} 

{population', position', vision', agentSugar', agentGulture', 

metabolism', children', diseases', agentimmunity', age', sex', initial Sugar') = (1) 
applyMating{rndNewSweep{potentialMatingPairs), population, position, vision, 
agentSugar, agentGulture, metabolism, children, 
diseases, agentimmunity, age, maxAge, sex, initialSugar) 


(1) Call applyMating function on the agents in sequence. 


4.11. Credit. 

_ M akeLoans _ 

A Agents 
EStep 

population' = population 

sex' = sex 

position' = position 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentGulture' = agentGulture 

agentimmunity' = agentimmunity 

diseases' = diseases 

children' = children 

metabolism' = metabolism 

initialSugar' = initialSugar 

{loanBook', agentSugar') = 

applyLoans{rndNewSweep{position), population, position, agentSugar, age, sex. 


(1) Call apply Loans on each agent in turn. 
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apply Loans : seq AGENT 
xFAGENT 

X AG ENT ^ POSITION 
xAGENT^N 
xAGENT^N 
xAGENT ^ SEX 
X {AGENT ^ {AGENT x (N x N))) 
xN 
-)■ 

{AGENT ^ {AGENT x (N x N))) 
xAGENT^N 


y population : ¥ AGENT; 
position : AGENT ^ POSITION; 
sex : AGENT ^ SEX; 
age : AGENT ^ N; 
agentSugar : AGENT -E- N; 
head, ag : AGENT; 

loanBook, loans : {AGENT ^ {AGENT x (M x M))); 
tail : seq AGENT; 
step : N; 

3 new AgentSugar : AGENT -H-M; newLoans : {AGENT -H- {AGENT x (N x N))); 
neighbours : P AGENT \ 

neighbours = {b : AGENT \ vonNeumanNeighbour{lender, b,position)} (1) 
{newLoans, new AgentSugar) = 

singleLender Loans{lender, as Seq {neighbours), 

agentSugar, age, sex, loanBook, step) 

applyLoans{{), population, position, agentSugar, age, sex, loanBook, step) = (2) 

{loanBook, agentSugar) 

apply Loans {{lender) ^ tail, population, position, 

agentSugar, age, sex, loanBook, step) = (3) 

apply Loans {tail, population, position, new AgentSugar, age, sex, newLoans, step) 


(1) Construct the set of neighbours of an agent lender, the updated loan book 
and the updated sugar levels gotten by the lender giving loans to its neigh¬ 
bours; 

(2) Base Case: Nothing to do just return existing values; 

(3) Recursive case: Recursively call applyLoans on the remainder of the agents 
(excluding the first agent lender) and the new loan and sugar levels gotten 
by lender generating new loans. 
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singleLenderLoans : AGENT 
X seq AGENT 
AGENT 
xAGENT^n 
xAGENT^n 
xAGENT ^ SEX 
X {AGENT FA {AGENT x (N x N))) 
xN 
-)■ 

{AGENT FA {AGENT x (N x N))) 
xAGENT^N 

^ sex : AGENT ^ SEX] 
age : AGENT ^ N; 
agentSugar : AGENT -A 
head, lender : AGENT ; 

loanBook : {AGENT FA {AGENT x (N x N))); 
tail : seq AGENT] 
step : N; 

singleLenderLoans{lender, (), agentSugar, age, sex, loanBook, step) = 

{loanBook, agentSugar) 

singleLenderLoans{lender, {head) ^ tail, agentSugar, age, sex, loans, step) = 
if canLend{age{lender), sex{lender), agentSugar {lender)) 

A willBorrow{age{head), sex{head), agentSugar {head), {head} <1 YUTL{loanBook))) 

then 

3 new AgentSugar : AGENT ha M; 

new Loans : {AGENT x {AGENT x (M x M))); amt : N | 

amt = min{amtAvail{age{lender), sex{lender), agentSugar {lender)), 
amtReq{agentSugar {head))) 

new AgentSugar = agentSugar © {lender i-A- agentSugar {lender) — amt, 
head i-A- agent Sugar {head) + amt} 

newLoans = loanBook U {{lender, {head, {amt, step + DURATION)} 
singleLender Loans{lender, tail, new AgentSugar, age, sex, newLoans, step) 

else 

singleLenderLoans{lender, tail, agentSugar, age, sex, loanBook, step) 

singleLender Loans calculates all loans that a particular agent can give to its 
neighbours. 

(1) If there are no loans in the sequence then just return the current loans and 
sugar levels as is; 

(2) If the loan sequence is not empty then apply the payment details to the 
hrst loan and make the payments on the rest: 
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a) If the first loan is capable of being paid by the borrower we simply move 
the correct amount of sugar from the borrower to the lender; 

b) If the borrower cannot pay off the loan then they pay back half their 
sugar and the loan is renegotiated for the remainder. 

Using these functions we can now specify the PayLoans part of the Credit rule. 

_ PayLoans _ 

AAgents 

EStep 

population' = population 

sex' = sex 

position' = position 

vision' = vision 

age' = age 

maxAge' = max Age 

agentCulture' = agentCulture 

agentimmunity' = agentimmunity 

children' = children 

diseases' = diseases 

metabolism' = metabolism 

initial Sugar' = initial Sugar 

3dueLoans,newLoans : F{AGENT x {AGENT x (N x N))) • 
dueLoans = loanBook > {iaxi{loanBook) > {a : (N x N) | second{a) = step}) (1) 
{newLoans, agentSugar') = makePayments{asSeq{dueLoans), 0, agentSugar){2) 
loanBook' = {loanBook \ dueLoans) U newLoans (3) 

(1) We create the set of due loans; 

(2) We now create the set of renegotiated loans and update the agentSugar 
levels using the makePayments function; 

(3) Finally we update the loan book by removing all loans that were due and 
adding any new renegotiated loans. 

5. Added Spice 

5.1. Introduction. We have dehned all the rules so far under the assumption 
that there is only one resource (known as sugar). The hnal rule, Trade, is only 
dehned for simulations with at least two resources. In fact the rules are meant to 
be general enough that they will work with any number of resources although we 
know of no sugarscape based simulation that used more than two resources. The 
second resource is known as spice. 

We will show how to extend the rules to deal with two resources. In order to 
avoid unnecessary clutter and make the differences as clear as possible we will 
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show the differences between the one and two resource schemas with boldface. 
Any part of a schema that is not an exact copy of the previously defined version 

will be in boldface. 

5.2. Basic Types. The basic types are copies of those already defined for sugar. 


MAXSPICEMETABOLISH : N (1) 

SPICEGROWTH ; N (2) 

MAX SPICE : N (3) 

INITIALSPICEMIN, INITIALSPICEMAX : M (4) 

SPICEPRODUCTION, SPIC EGON SUMPTION : N(5) 
SPICECOMB AT LIMIT : M ( 6 ) 

SPICECHILDAMT : M (7) 


(1) Agents metabolise spice during each move at an individually set rate less 
than MAXSPICEMETABOLISH] 

(2) Spice grows back at a predefined rate; 

(3) Each location can hold a set maximum amount of spice; 

(4) Agents created after mating start with an initial spice endowment; 

(5) Pollution can be caused by production and consumption of spice; 

( 6 ) SPICECOMBATLIMIT is required to help determine the reward from 
attacking an agent using the combat rule. 

(7) We posit that a minimum amount of spice is needed for agent mating to 
occur. 

Note that these constants are replicas of their sugar counterparts. 

5.3. The SpiceScape. The spice grid contains everything in the Lattice scheme 
and just adds information on the extra spice resource. 

_ SpiceLattice _ 

Lattice 

spice : POSITION ^ N 
maxSpice : POSITION N 

dom spice = dom maxSpice = POSITION (1) 

Vx : POSITION • spice(x) < maxSpice(x) < MAXSPICE (2) 


(1) Every location has an associated amount of spice and maximum carrying 
capacity; 

(2) Every position’s spice levels are within the acceptable levels. 


5.4. Agents. 
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_ SpiceAgents _ 

Agents 

agentSpice : AGENT -E- N 
initialSpice : AGENT -E- N 
spiceMetabolism ; AGENT -i^ N 
spiceLoanBook : AGENT (AGENT x (N x N)) 

doni spiceMetabolism = doni agentSpice == dominitialtSpicepopulation(l) 
doni spiceLoanBook C population (2) 

doni(ran spiceLoanBook) C population 
Vx : AGENT • x G population ^ 

spiceMetabolism(x) < MAXSPICEMETABOLISM (3) 


(1) Every agent has a spice metabolism and a spice store; 

(2) The spiceLoanBook has the same invariants as the original loanBook; 

(3) Every agents metabolism is less than or equal to the dehned maximum. 


Finally we combine them into an overall schema as before: 


_ Spice Scape _ 

SpiceAgents 

SpiceLattice 

Step 


The initialisation scheme and tick schemas are also largely unchanged. 
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_ Initial Spice Scape _ 

SpiceAgents' 

SpiceLattice' 

Step' 

step' = 0 

# population' = INITIALPOPULATIONSIZE 
loanBook' = 0 
spiceLoanBook' = 0 
Va: AGENT • 

a G population' ^ 
age'{a) = 0 
diseases' (a) = 0 
children' {a) = 0 

INITIALSUGARMIN < agent Sugar'{a) < INITIALSUGARMAX 
initial Sugar'(a) = agentSugar' {a) 

INITIALSPICEMIN < agentSpice'(a) < INITIALSPICEMAX 
initialSpice'(a) = agentSpice'(a) 
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L spice - 

ASpiceAgents 
A Step 

population' = population 

position' = position 

sex' = sex 

vision' = vision 

maxAge' = maxAge 

metabolism' = metabolism 

initial Sugar' = initial Sugar 

initialSpice' = initialSpice 

spiceMetabolism' = spiceMetabolism 

agentCulture' = agentCulture 

children' = children 

loanBook' = loanBook 

spiceLoanBook' = spiceLoanBook 

agent Immunity' = agentimmunity 

diseases' = diseases 

step' = step + 1 

Vx : AGENT • x G population ^ 

{age'{x) = age{x) + 1 

A agentSugar'{x) = agentSugar{x) — metabolism{x) 

A agentSpice'(x) = agentSpice(x) — spiceMetabolism(x)) 


5.5. Rules. As well as defining the final rule, Trade, we will also expand the other 
rules to allow them to operate on a simulation with two resources. We dehne the 
new rule {Trade) first. 

5.6. Agent Trade T. 

Agent Trade T: 

• Agent and neighbour compute their MRSs; if these are equal then end, 
else continue; 

• The direction of exchange is as follows: spice flows from the agent with 
the higher MRS to the agent with the lower MRS while sugar goes in 
the opposite direction; 

• The geometric mean of the two MRSs is calculated-this will serve as 
the bargaining price, p; 

• The quantities to be exchanged are as follows: if p>l the p units of 
spice for 1 unit of sugar; if p < 1 the 1 /p units of sugar for 1 unit of 
spice; 
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• If this trade will (a) make both agents better off (increases the welfare 
of both agents), and (b) not cause the agents’ MRSs to cross over one 
another, then the trade is made and return to start, else end. 


MRS is calculated simply for an agent as the fraction obtained by dividing its 
spice level times its sugar metabolism by its spice metabolism times its sugar level, 
as set out below. 


MRS 

V sugar, sugar Metabolism, spice, spiceM etabolism : N • 
MRS {sugar, sugar Metabolism, spice, spiceM etabolism) = 
{spice * sugar Metabolism) / {spiceM etabolism * sugar) 


tradePairs constructs a sequence of all possible trading partners based on the 
proximity of the agents to each other. 


tradePairs : seq AGENT x AGENT POSITION ^ seq{AGENT x AGENT) 

ytail : seq AG ENT; positions : AGENT POSITION; a : AGENT • 
tradePairs{{), positions) = () 
tradePairs{{a) ^ tail,positions) = 

trade{tail,positions asSeq{{b ; agent \ adjacent{position{a),position{b)) • {a,b)}) 
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_ T rade _ 

ASpiceAgents 


spiceLoanBook' = spiceLoanBook 
loanBook' = loanBook 
population' = population 
initial Sugar' = initial Sugar 

initialSpice' = initialSpice 

sex' = sex 

metabolism' = metabolism 
spiceM etabolism' = spiceM etabolism 
position' = position 


vision' = vision 


age' = age 
maxAge' = maxAge 
agentCulture' = agentCulture 
agentimmunity' = agentimmunity 
children' = children 

3allPairs : seq{AGENT, AGENT) \ allPairs = tradePairs{asSeq{population), position) 
{agentSugar', agentSpice') = (1) 

executeTrades{allTrades{chooseExclusiveTrades{allPairs), 

{agentSugar, agentSpice), metabolism, spiceM etabolism)) 


Trade is similar to Mating in that trading must be done in exclusive pairs. An 
agent cannot carry out two simultaneous trades and the rule forces each agent 
to trade with all its neighbours in some sequence. As with Mating we construct 
conflict free sets of trading pairs that can proceed simultaneously and then order 
these sets. 


chooseExclusiveTrades : AGENT x AGENT ^ seqF{AGENT, AGENT) 

\/a,b: AGENT] trading Pairs : AGENT ^ AGENT • 
chooseExclusiveTrades{0) = () 
chooseExclusiveT rades {trading Pairs) = 

3 maxSet : AGENT -H- AGENT \ maxSet C tradingPairs 
A ((a, b) G tradingPairs A (a, b) ^ maxSet) 

3 c : AGENT \ {(a, c), (c, a), (6, c), (c, &)} H maxSet ^ 0) 

{maxSet) ^ chooseExclusiveTrades{tradingPairs \ maxSet) 
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executeTrades : seq{AGENT x AGENT) x AGENT ^Nx AGENT ^ N 
^ {AGENT ^ N, AGENT ^ N) 

ytail : seq{AGENT x AGENT); head : {AGENT x AGENT); 
agentSugar, agentSpice, metabolism, spiceMetabolism : AGENT N • 
executeTrades{{), agentSugar, agentSpice, metabolism, spiceMetabolism) = 
{agentSugar, agentSpice) 

executeTrades {{head) tail, agentSugar, sugar) = 

3 new AgentSugar, new AgentSpice : AGENT -H- M | {new AgentSugar, new AgentSpice) 
allTrades{head, {agentSugar, agentSpice, metabolism, spiceM etabolism) 

• executeTrades{tail, new AgentSugar, new Agent Spice, metabolism, spiceM etabolism) 


allTrades recursively goes through the sequence of trading partners and gets 
each trading pair to update the sugar and spice levels based on their trades. 


allTrades : seq{AGENT x AGENT) x {{AGENT N) x {AGENT N))x 
{AGENT N) X {AGENT N) 

^{{AGENT N) X {AGENT N)) 

yhead : AGENT x AGENT; tail : seq{AGENT x AGENT) 
sugar, spice, sugar Metabolism, spiceM etabolism : AGENT h-> N • 
allTrades{{), {sugar, spice), sugarMetabolism, spiceMetabolism) = {sugar, spice) 
allTrades {{head) ^ tail, {sugar, spice), sugar Metabolism, spiceMetabolism) = 

allTrades{tail, pairTrade{head, {sugar, spice), sugarMetabolism, spiceM etabolism) 
, sugarMetabolism, spiceMetabolism) 


Each trading partnership will execute a series of trades until their MRS scores 
cross over. pairTrade is complicated by the fact that there are multiple options: 

(1) If their MRS scores are equal then they perform no trades; 

(2) If their MRS scores are not equal then the direction of trade will depend 
on which MRS score is higher; 

a) Within a trade a value of p (based on MRS scores) determines the price 
of the resources. 

(3) Trades between the pair continue until their new and old MRS scores cross 


over. 
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pairTrade : {AGENT x AGENT) x {{AGENT M) x {AGENT M))x 
{AGENT N) X {AGENT N) 

^{{AGENT N) X {AGENT N)) 

a,b AGENT] sugar, spice : AGENT >+^ N; 
metabolism, spiceMetabolism ; AGENT N | 

3p,mrsA,mrsB,newMrsA,newMrsB : A; new Sugar, new Spice : AGENT h-> N 
mrsA = MRS{sugar{a), metabolism{a), spice{a), spiceMetabolism{a)) 
mrsB = MRS{sugar{b),metabolism{b), spice{b), spiceMetabolism{b)) 
p = y/mrsA * mrsB 
{mrsA > mrsB A p > 1) ^ 

{newSugar = sugar © {(a, sugar{a) + 1), {b, sugar{b) — 1)} 

A newSpice = spice © {a i—)■ spice{a) — p,b\-A- spice{b) + p}) 

{mrsA > mrsB A p < 1) ^ 

{newSugar = sugar © {a i—)■ sugar{a) + (1 div p),b\-A sugar{b) — (1 div p)} 
newSpice = spice © {a i-A- spice{a) — 1, 61 —)■ spice{b) + 1}) 

{mrsA < mrsB A p > 1) ^ 

{newSugar = sugar ©{fen-)- sugar{b) + 1, a i-A- sugar{a) — 1} 

A newSpice = spice © {6 i-A- spice{b) — p,ae^ spice{a) + p}) 

{mrsA < mrsB A p < 1) ^ 

{newSugar = sugar © {6 h-)■ sugar {b) + (1 div p),ae^ sugar {a) — (1 div p)} 

A newSpice = spice © {6 i-A- spice{b) — 1, a i—)■ spice{a) + 1}) 
new Mrs A = MRS{newSugar{a),metabolism{a),newSpice{a), spiceM etabolism{a)) 
new MrsB = M RS {new Sugar {b), metabolism{b), new Spice{b), spiceM etabolism{b)) 

• 

pairTrade{{a,b), {sugar, spice), metabolism, spiceMetabolism) = {sugar, spice) 

AA mrsA = mrsB 

pairTrade{{a,b), {sugar, spice), metabolism, spiceMetabolism) = {new Sugar, new Spice) 
AA {{mrsA > mrsB A newMrsA < newMrsB) 

V {mrsA < mrsB A newMrsA > newMrsB)) 
pairTrade{{a, b), {sugar, spice), metabolism, spiceM etabolism) = 

pairTrade{{a, b), {newSugar, new Spice), metabolism, spiceM etabolism) 

AA {{mrsA > mrsB A newMrsA > newMrsB) 

V {mrsA < mrsB A newMrsA < newMrsB)) 


5.7. Asynchronous Trade. 
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_ Trade _ 

ASpiceAgents 

spiceLoanBook' = spiceLoanBook 
loanBook' = loanBook 
population' = population 
sex' = sex 

metabolism' = metabolism 
spiceM etabolism' = spiceMetabolism 
initialSugar' = initialSugar 
initialSpice' = initialSpice 
position' = position 
vision' = vision 
age' = age 
maxAge' = maxAge 
agentCulture' = agentCulture 
a gent Immunity' = agentimmunity 
children' = children 
3 traders ; F{AGENT x AGENT) • 

{agentSugar',agentSpice') = (1) 

allTrades(tradePairs{rndNewSweep{position), positions), {agentSugar, agentSpice), 
metabolism, spiceM etabolism) 


(1) The new sugar and spice allocations are derived by conducting all possible 
trades using the recursive helper function allTrades. 

5.7.1. Growback. 

_ Growback spice _ 

ASpiceLattice (1) 

pollution' = pollution 
maxSugar' = maxSugar 

maxSpice' = maxSpice (2) 

sugar' = sugar © {x : POSITION • 

X min{{sugar{x) + SUGARGROWTH, maxSugar{x)})} 

spice' = spice © {x : POSITION • 

X I—)■ min({spice(x) + SPICEGROWTH, maxSpice(x)})} (3) 


(1) Lattice is replaced with SpiceLattice. In all subsequent schemas SpiceLattice 
replaces Lattice and SpiceAgent replaces Agent, 

(2) maxSpice remains unchanged; 

(3) The new spice levels are calculated using the same simple formula used for 
sugar growback. 
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5.7.2. Seasonal Growback. 

_ SeasonalGrowhack spice _ 

ASpiceLattice 

EStep 

pollution' = pollution 
maxSugar' = maxSugar 

maxSpice' = maxSpice 

Vx : POSITION • 

{step div SEASONLENGTH) mod 2 = 0^ sugar' = 

{x : POSITION I first{x) < M div 2 • 

X min{{sugar{x) + SUGARGROWTH, max Sugar {x)})} 

U 

{x ; POSITION I first{x) > M div 2 • 

X 1 -^ min{{sugar{x) + SUGARGROWTH div WINTERRATE, maxSugar{x)})} 
Vx : POSITION • 

{step div SEASONLENGTH) mod 2^0^ sugar' = 

{x : POSITION I first{x) < M div 2 • 

X ^ min{{sugar{x) + SUGARGROWTH div WINTERRATE, max Sugar {x)})} 
U 

{x : POSITION-, y:n\ first{x) > M div 2 • 

X 1-4 min{{sugar{x) + SUGARGROWTH, max Sugar {x)})} 

Vx ; POSITION • (1) 

(step div SEASONLENGTH) mod 2 = 0^ spice' = 

{x ; POSITION I first(x) < M div 2 • 

X 1-4 min({spice(x) + SPICEGROWTH, maxSpice(x)})} 

U 

{x : POSITION I first(x) > M div 2 • 

X ^ min({spice(x) + SPICEGROWTH div WINTERRATE, maxSpice(x)})} 
Vx : POSITION • 

(step div SEASONLENGTH) mod 2^0^ spice' = 

{x : POSITION I first(x) < M div 2 • 

X ^ min({spice(x) + SPICEGROWTH div WINTERRATE, maxSpice(x)})} 

U 

{x ; POSITION I first(x) > M div 2 • 

X 1-4 min({spice(x) + SPICEGROWTH, maxSpice(x)})} 


(1) Seasonal growback adds a rule for spice grow back that is an exact replica 
of the sugar rule. We note that we only use the one WINTERRATE 
instead of a separate rate for sugar and spice. In the absence of any explicit 
direction on this point this solution seems to be the most obvious. 
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5.8. Movement. In order to update Movement we will need to implement a wel¬ 
fare function that can be used to measure the desirability of a location. With two 
resources the desirability of any location becomes a subjective measure, what one 
agent may rate highly another may not. Welfare is dependent on the agents current 
levels of spice and sugar, so an agent with low spice levels may consider a location 
containing spice more desirable than one containing sugar. Overall the desirability 
of a location is determined by the agents current resource levels (wealth) and the 
relative metabolism rates for each resource. 

This is in contrast to the previous approach where welfare just equaled the 
amount of sugar in a location. This welfare measure is precisely dehned in the 
book and we follow this dehnition. 


welfare ;NxNxMxMxMxN^A 

V agentSugar, sugar Metabolism, 
agentSpice, spiceMetabolism, 
locationSugar, locationSpice : N • 

welfare{agentSugar, sugar M etabolism, agentSpice, spiceM etabolism, 
locationSugar, locationSpice) = 

{locationSugar + agentSugar) * {sugarMetabolism div 
{sugar Metabolism + spiceM etabolism)) 

^{locationSpice + agentSpice) * {spiceMetabolismdiv 
{sugar Metabolism -|- spiceM etabolism)) 


Movement can now be restated by replacing the previous measure of a locations 
desirability (sugar level) with this new measure and the updating of spice levels. 
In all other respects the schema remains unchanged. 
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— A4 ov €-TncTit]jQ^gi(^gpi(^^ - 

ASpiceScape 

step’ = step 

population' = population 
maxSugar' = maxSugar 
pollution' = pollution 
sex' = sex 
vision' = vision 
age' = age 

initial Sugar' = initial Sugar 

initialSpice' = initialSpice 

metabolism' = metabolism 

maxAge' = maxAge 

agentCulture' = agentCulture 

loanBook' = loanBook 

diseases' = diseases 

agent Immunity' = agentimmunity 

children' = children 

maxSpice' = maxSpice 

spiceLoanBook' = spiceLoanBook 

spiceMetabolism' = spiceMetabolism 

Va : AGENT- I ; POSITION • 

a G population ^ distance{position'{a),position{a)) < vision{a) 
distance{position{a), 1) < vision{a) A (/ 0 ran position')) 

welfare(agentSugar(a), metabolism(a), (1) 

agentSpice(a), spiceMetabolism(a), sugar(l), spice(l)) 

< 

welfare(agentSugar(a), metabolism(a), 
agentSpice(a), spiceMetabolism(a), 

sugar (position' (a)), spice(position' (a))) 

agentSugar' = {a ; AGENT \ a G population • 
a i-G- agentSugar{a) + sugar {position' (a))} 
sugar' = sugar © {loc : POSITION \ loc G lanposition' • loc i-G- 0} 

agentSpice' = {a : AGENT | a G population • (2) 

a I—)■ agentSpice(a) 

©spice (position' (a))} 

spice' = spice © {loc ; POSITION | loc G ran position' • loc i-© 0} (3) 

(1) This is a copy of the original proposition with the welfare function now 
replacing the previous sugar level check; 

(2) Agents consume spice at their new locations; 

(3) Locations with agents present now have no remaining spice. 
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5.9. Pollution Formation Pn,x- The M ovementspicePoiiuUon schema has the same 
alterations as the MovementhasicSpice 
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— M ovcTficntpQii^iiQjigpif,^ - 

ASpiceScape 
step' = step 

population' = population 
maxSugar' = maxSugar 
sex' = sex 

pollution' = pollution 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentCulture' = agentCulture 

loanBook' = loanBook 

children' = children 

agent Immunity' = agent Immunity 

initial Sugar' = initial Sugar 

initialSpice' = initialSpice 

diseases' = diseases 

metabolism' = metabolism 

spiceMetabolism' = spiceMetabolism 

spiceLoanBook' = spiceLoanBook 

maxSpice' = maxSpice 

Va : AGENT- I : POSITION \ a G dom{position') • 

a G population ^ distance{position'{a),position{a)) < vision{a) 
{distance{position{a), 1) < vision{a) A (/ ^ ran position')) 

welfare(agentSugar(a), metabolism(a), 
agentSpice(a), spiceMetabolism(a), sugar(l), spice(l)) 
div(l + pollution(l)) 

< 

welfare(agentSugar(a), metabolism(a), 

agentSpice(a), spiceMetabolism(a), sugar (position'(a)), 
spice(position'(a))) div (1 + pollution(position'(a))) 

agentSugar' = {a : AGENT \ a G population • 
a i-G- agentSugar{a) + sugar {position' (a))} 
sugar' = sugar 0 {/ : POSITION \ I G ranposition' • / i-G- 0} 

agentSpice' = {a : AGENT | a G population • 
at—)■ agentSpice(a) + spice(position'(a))} 
spice' = spice © {1 ; POSITION | 1 G ran position' • 1 i-G- 0} 
pollution' = pollution® 

{I : POSITION-, X : AGENT \ position'{x) = I • 

I ^ {PRODUGTION * sugar{l) + GONSUMPTION * metabolism{x))} 
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5.10. Pollution Diffusion. 

_ PollutionDif fusion spice _ 

ASpiceLattice 

EStep 

maxSugar' = maxSugar 
sugar' = sugar 

maxSpice' = maxSpice 
spice' = spice 

{step mod POLLUTION RATE 7 ^ 0) pollution' = pollution 

{step mod POLLUTION RATE = 0) pollution' = 

{I : POSITION • I i-A {pollution{north{l)) + pollution{south{l)) 
+pollution{east{l)) + pollution{west{l))) div 4} 


5.11. Replacement. 


ASpiceAgents 

population' = population \ {a : AGENT \ age{a) = maxAge{a) 

V agent Sugar {a) = OV agentSpice(a) = 0 • a} 
loanBook' = population' < loanBookt> 

{x : AGENT x (N x N) | first{x) G population'} 

spiceLoanBook' = population' < spiceLoanBook> 

{x ; AGENT x (N x N) | first(x) G population'} 

Va: AGENT • 

a G population' ^ 

{sex{a) = sex'{a) A vision{a) = vision'{a) 

A maxAge{a) = maxAge'{a) A agentGulture{a) = agentGulture'{a) 
A position{a) = position'{a) A age{a) = age'{a) 

A agentSugar{a) = agentSugar'{a) 

A metabolism'{a) = metabolism{a) 

A diseases' {a) = diseases{a) 

A agentimmunity'{a) = agentlmmunity{a) 

A children'{a) = children{a) 

A initial Sugar'{a) = initial Sugar {a) 

A agentSpice'(a) = agentSpice(a) 

A initialSpice'(a) = initialSpice(a) 

A spiceMetabolism'(a) = spiceMetabolism(a)) 
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_ Replacement spice _ 

ASpiceAgents 

# population' = INITIALPOPULATIONSIZE 
loanBook' = loanBook 
spiceLoanBook' = spiceLoanBook 
Va :AGENT • 
a G {population) 

{a G population' 

A sex{a) = sex'{a) A vision{a) = vision'{a) 

A maxAge{a) = maxAge'(a) A agentCulture{a) = agentCulture'{a) 

A position{a) = position'{a) A age{a) = age'{a) 

A agentSugar'{a) = agentSugar{a) 

A metabolism'{a) = metabolism{a) 

A initial Sugar'{a) = initial Sugar {a) 

A initialSpice'(a) = initialSpice(a) 

A agentSpice'(a) = agentSpice(a) 

A spiceMetabolism'(a) = spiceMetabolism(a) 

A diseases' {a) = diseases{a) 

A agent Immunity' {a) = agent Immunity {a) 

A children'{a) = children{a)) 

Ma:AGENT • 

a G population' \ population ^ {age'{a) = 0 

A INITIALSUGARMIN < agentSugar'{a) < IN IT I ALSU G ARM AX 
A initial Sugar'{a) = agentSugar'{a) 

A INITIALSPICEMIN < agentSpice'(a) < INITIALSPICEMAX 
A initialSpice'(a) = agentSpice'(a) 

A diseases' {a) = 0 A children' {a) = 0 ) 


5.12. Agent Mating. 
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_ Agent M ating _ 

ELattice 
A Agents 

loanBook' = loanBook A spiceLoanBook' = spiceLoanBook 

3 potential Mating Pairs : F{AGENT x AGENT) \ (1) 

potential Mating Pairs = {(a ; AGENT,h : AGENT) \ sex{a) ^ sex{h) 

A isEertile{age{a), sex{a)) A isEertile{age{head), sex{head)) 

A adjacent{position{a),position{head))} 

{population', position', vision', agentSugar', agentGulture', metabolism' 

, children', diseases', agentimmunity', age', sex', initial Sugar', 
spiceMetabolism', agentSpice', initialSpice') = (2) 

concurrentMating{getGon fictEreePairs{potentialMatingPairs), population, position, vision, 
agentSugar, agentGulture, metabolism, children, 
diseases, agentimmunity, age, maxAge, sex, initialSugar, 
spiceMetabolism, agentSpice, initialSpice) 


(1) Generate the set of all possible mating pairs; 

(2) Recursively proceed with concurrent mating within the conflict free subsets. 
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concurrent Mating : sec[¥{AGENT x AGENT) 
AGENT 

X AG ENT ^ POSITION 
xAGENT 
xAGENT^N 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT ^ P AGENT 
X AG ENT -E-Pseq-B/T 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT 
xAGENT ^ SEX 
xAGENT^N 
X AGENT 
X AGENT 
X AGENT 

PAGE AT 

X AG ENT ^ POSITION 
xAGENT 
xAGENT^N 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT ^ P AGENT 
X AG ENT -H-PseqEJT 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT 
xAGENT ^ SEX 
xAGENT^n 
X AGENT 
X AGENT 
X AGENT 

\/tail : seqP(AGEAT x AGENT); 

head : P(AGEAT x AGENT); 

population : P AGENT; 

position : AGENT ^ POSITION; 

vision : AGENT -H- Ni; 

agentSugar, agentSpice ; AGENT -H- M; 

agentGulture ; AGENT -+)■ seqBIT; 

metabolism, spiceMetabolism : AGENT -H- M; 

children : AGENT ^ PAGE AT; 

diseases : AGENT -H-PseqE/T; 

agentimmunity : AGENT -H-seqE/T; 

age : AGENT ^ N; 

maxAge : AGENT Ni; 

sex : AGENT ^ SEX; 

initialSugar, initialSpice ; AGENT N; 
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3newpopulation : F AGENT] 
newposition : AGENT ^ POSITION] 
newvision : AGENT Ni; 

newagentSugar, newagentSpice : AGENT N; 

newagentGulture : AGENT -H- seqBIT] 

newmetabolism, newspiceMetabolism ; AGENT -H- M; 

newchildren : AGENT F AGENT] 

newdiseases : AGENT -H-Pseq-BIT; 

newagent Immunity : AGENT -^seqBIT] 

new age ; AGENT -H- N; 

newmaxAge : AGENT Ni; 

newsex : AGENT -H- SEX] 

newinitialSugar, newinitialSpice : AGENT N; | 

{newpopulation, newposition, newvision, newagentSugar, newagentGulture, 

newmetabolism, newchildren, newdiseases, new agent Immunity, newage, newmaxAge, 
newsex, newinitial Sugar, newspiceMetabolism, newagentSpice, newinitialspice) = 
applyMating{asSeq{head), population, position, vision, 
agentSugar, agentGulture, metabolism, children, 
diseases, agent Immunity, age, maxAge, sex, initialSugar, 
spiceMetabolism, agentSpice, initialspice) • 
concurrentMating{{), population, position, vision, 
agentSugar, agentGulture, metabolism, children, 
diseases, agent Immunity, age, maxAge, sex, initialSugar) = 

{population, position, vision, agentSugar, agentGulture, metabolism, 
children, diseases, agentimmunity, age, maxAge, sex, initialSugar, 
spiceMetabolism, agentSpice, initialspice) 

concurrentMating{{head) ^ tail, population, position, vision, 
agentSugar, agentGulture, metabolism, children, 
diseases, agentimmunity, age, maxAge, sex, initialSugar) = 
concurrentM ating {tail, newpopulation, newposition, newvision, 

newagentSugar, newagentGulture, newmetabolism, newchildren, 
newdiseases, new agent Immunity, newage, newmaxAge, newsex, newinitial Sugar, 
newspiceMetabolism, newagentSpice, newinitialspice) 
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apply Mating : seq{AGENT x AGENT) 
AGENT 

X AG ENT ^ POSITION 
xAGENT 
xAGENT^N 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT ^ P AGENT 
X AG ENT -E-Pseq-B/T 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT 
xAGENT ^ SEX 
xAGENT^N 
X AGENT 
X AGENT 
X AGENT 

PAGE AT 

X AG ENT ^ POSITION 
xAGENT 
xAGENT^N 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT P AGENT 
X AG ENT -H-PseqEJT 
X AGENT ^seqBIT 
xAGENT^N 
X AGENT 
xAGENT ^ SEX 
xAGENT^N 
X AGENT 
X AGENT 
X AGENT 


\/ population : PAGE AT; 

position : AGENT ^ POSITION- 

sex : AGENT SEX-, 

vision -. AGENT -H- Ni; 

age : AGENT N; 

initialSugar, initialSpice ; AGENT N; 
maxAge : AGENT -H- Ni; 

metabolism, spiceMetabolism : AGENT N; 
agentSugar, agentSpice : AGENT -H- N; 
agentGulture : AGENT seqBIT-, 
children : AGENT ^ PAGE AT; 
agentimmunity : AGENT -H-seqE/T; 
diseases : AGENT -H-PseqE/T; 
head : AGENT x AGENT-, 
tail -. seq(AGEAT x AGENT)-, • 
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3 offspring, a, b : AGENT ; 
newsex : AGENT -H SEX] 
newvision : AGENT Ni; 

newmetabolism, newagentSugar, newinitialSugar : AGENT M; 
newspiceMetabolism, newagentSpice, newinitialSpice ; AGENT N; 

newmaxAge : AGENT Ni; 
newagentGulture : AGENT -H- seqBIT; 
newchildren : AGENT ¥ AGENT] 
newagentImmunity : AGENT -Hseq-B/T; 
inkeritedimmunity : seqBIT] 
inkeritedGulture : seqBIT] 

I offspring ^ population 
a = firstijiead) A b = second{kead) 

newckildren : ckildren U {of f spring i—)■ 0 a i—)■ ckildren{a) U [of f spring}, 
b i-A ckildrenib) U {of f spring}} 
newsex = sex U {of f spring i—)■ male} 

V newsex = sex U {offspring i—)■ female} 
newvision = vision U {offspring i-A vision{a)} 

V newvision = vision U {of f spring i—)■ vision{b)} 
newmaxAge = max Age U {of f spring i-A maxAge{a)} 

V newmaxAge = max Age U {of f spring i-A maxAge{b)} 
newmetabolism = metabolism U {of f spring i—)■ metabolism{a)} 

V newmetabolism = metabolism U {of f spring i-A metabolism{b)} 
newspiceMetabolism = spiceMetabolism U {offspring i—)■ spiceMetabolism(a)} 

V newspiceMetabolism = spiceMetabolism U (offspring i-A spiceMetabolism(b)} 
newinitialSugar = initialSugar® 

{of f spring i-A initial Sugar {a)/2 + initialSugar{b)/2} 
newagentSugar = agentSugar © {of f spring i-A initialSugar{of f spring), a i-A initialSugar{a)/2,1 
newinitialSpice = initialSpice© 

(offspring i-A initialSpice(a)/2 + initialSpice(b)/2} 
newagentSpice = agentSpice U (offspring i—)■ initialSpice(offspring), a ha initialSpice(a)/2, f 
A Vn : 1 .. IMMUNITYLENGTH • 

{inkeritedlmmunity{n) = agentlmmunity{a){n) 

V inkeritedlmmunity{n) = agentlmmunity{b){n)) 
newagentimmunity : agentimmunity U {of f spring ha inkeritedimmunity} 

A Vn : 1 .. GULTUREGOUNT • 

{inker itedGulture{n) = agentGulture{a){n) 

V inkeritedGulture{n) = agentGulture{b){n)) 
newagentGulture : agentGulture U {of f spring ha inker itedGulture} 
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applyMating{{), population, position, vision, agentSugar, agentCulture, 

metabolism, children, diseases, agentimmunity, age, maxAge, sex, initialSugar, 

spiceMetabolism, agentSpice, initialSpice) = 

{population, position, vision, agentSugar, agentCulture, metabolism, 
children, diseases, agentimmunity, age, maxAge, sex, initialSugar, 

spiceMetabolism, agentSpice, initialSpice) 

apply Mating {{head) ^ tail, population, position, vision, agentSugar, agentCulture, 
metabolism, children, diseases, agentimmunity, age, maxAge, sex, initialSugar, 
spiceMetabolism, agentSpice, initialSpice) = 
if((3/oc : POSITION \ {adjacent{loc,position{ag))) 

V adjacent{loc,position{head)) A loc^ domposition) 

A (agentSpice(head) > initialSpice(head)) A (agentSpice(ag) > initialSpice) 

apply Mating {tail, population U {of f spring}, position U {of f spring i—)■ loc}, 
newvision, newagentSugar, newagentCulture, new metabolism, 
newchildren, diseases U {of f spring i—)■ 0}, new agent Immunity, 

age U {of f spring i—)■ 0}, newmaxAge, newsex, initialSugar, newspiceMetabolisn 

else 

apply Mating {tail, population, position, vision, agentSugar, agentCulture, 

metabolism, children, diseases, agentimmunity, age, maxAge, sex, initialSugar, 

spiceMetabolism, agentSpice, initialSpice) 


5.13. Culture. 
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— tlltUT spice - 

ASpiceAgents 

population' = population 

sex' = sex 

position' = position 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentSugar' = agentSugar 

children' = children 

loanBook' = loanBook 

diseases' = diseases 

initialSugar' = initialSugar 

metabolism' = metabolism 

a gent Immunity' = agentimmunity 

agentSpice' = agentSpice 

spiceMetabolism' = spiceMetabolism 

spiceLoanBook' = spiceLoanBook 

initialSpice' = initialSpice 

V a : AGENT • a G population ^ 

agentCulture'{a) = flipTags{agentCulture{a), 

asSeq{{b : AGENT \ adjacent{position{a),position{b))}), agentGulture) 


5.14. Disease. 
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_ ImmuneResponsespice _ 

ESpiceLattice 

ASpiceAgents 

EStep 

loanBook' = loanBook 
spiceLoanBook' = spiceLoanBook 

population' = population 

sex' = sex 

position' = position 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentCulture' = agentCulture 

diseases' = diseases 

children' = children 

agentSugar' = agentSugar 

initial Sugar' = initial Sugar 

initialSpice' = initialSpice 

metabolism' = metabolism 

agentSpice' = agentSpice 

spiceMetabolism' = spiceMetabolism 

agentimmunity' = {a : AGENT \ a G population • 

a i-G- applyDiseases{agentlmmunity{a), asSeq{diseases(a)))} 

'i X : AGENT • X E population ^ agentSugar'{x) = agentSugar {x) — 
#{(i : seqBIT \ d G diseases{a) A subseq{d, agentlmmunity{a)} 
'iX : AGENT • x E population agentSpice'{x) = agentSpice{x) — 
#{(i : seqBIT \ d G diseases{a) A subseq{d, agentlmmunity{a)} 
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_ Transmission spice _ 

ASpiceAgents 

loanBook' = loanBook 

initialSugar' = initialSugar 

spiceLoanBook' = spiceLoanBook 

population' = population 

sex' = sex 

position' = position 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentCulture' = agentCulture 

agentimmunity' = agentimmunity 

children' = children 

agentSugar' = agentSugar 

agentSpice' = agentSpice 

spiceMetabolism' = spiceMetabolism 

initialSpice' = initialSpice 

metabolism' = metabolism 

V a : AGENT • a G population 

diseases' {a) = diseases{a)U 

newDiseases{asSeq{visibleAgents{a,position, 1)), diseases) 


5.15. Inheritance. 
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_ Inheritance spice _ 

ASpiceAgents 

population' = population A sex' = sex 

position' = position A vision' = vision 

age' = age A maxAge' = maxAge 

agentCulture' = agentCulture A children' = children 

metabolism' = metabolism A spiceMetabolism' = spiceMetabolism 

diseases' = diseases A agentimmunity' = agentimmunity 

initialSugar' = initialSugar 

initialSpice' = initialSpice 

3 dying : F AGENT] 

inheritFromFemale,mheritFromMale : AGENT -E- (M x N) • (1) 

dominheritFromFemale = dominheritFromMale = population \ dying 
dying = {x : AGENT \ x G population A {age{x) = maxAge{x) V 

agentSugar{x) = OV agentSpice(x) = 0)} (2) 

V X : AGENT ; n,m : N \ x ^ population \ dying • 
getMother {x, children, sex) ^ dying ^ 
inheritFromFemale{x) = (0, 0) 
getFather{x, children, sex) ^ dying ^ 
inheritFromMale{x) = (0, 0) 

3m : AGENT \ m = getMother{x, children, sex) A m G dying 
inheritFromFemale{x) = 

{agentSugar{m) div ^{population n children{m) \ dying)) 
(agentSpice(m) div 7 ^(population fl children(m) \ dying)) (3) 

3 d : AGENT \ d = getFather{x, children, sex) A d E dying ^ 
inheritFromM ale{x) = 

agentSugar{d) div ^{population fl children{d) \ dying)) 

(agentSpice(m) div 7 )^(population fl children(ni) \ dying)) 

X G dying 

{agentSugar'{x) = OA agentSpice'(x) = 0) 

X ^ dying 

^ {agentSugar'{x) = agentSugar{x) 

+first{inheritFromMale{x)) + first{inheritFromFemale{x)) 

A agentSpice'(x) = agentSpice(x) 

+second(inheritProniFeniale(x)) + second(inheritProniMale(x))) 

loanBook' = disperseLoans{loanBook, asSeq{dying), children) 

spiceLoanBook' = 

disperseLoans(spiceLoanBook, asSeq(dying), children) (4) 


(1) Agents can now inherit two amonnts, a sugar inheritance and a spice 
inheritance; 
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(2) Death now occurs if either resource reaches zero; 

(3) The individual spice inheritance is calculated in the same way as the sugar 
inheritance; 

(4) Spice loans are dispersed amongst children. 


5.16. Combat. Combat is dehned only in terms of sugar. We can either accept 
this and assume combat is based only on sugar levels or we can extend combat 
by dehning new versions of wealth and reward. We note that no simulations 
combining combat with more than one resource are presented in the book. 

We can extend the combat rule with a few simple assumptions. First the wealth 
of an agent is used to determine if we can attack that agent or if an agent can 
retaliate against us. In the single resource scenario we simply used the sugar that 
an agent carried. With two resources we need to combine both sugar and spice. 
The simplest approach is to add these two amounts together and in the absence 
of any guidelines this seems the sensible option. 

availMoves requires only minor changes to return the set of all safe moves that 
an agent can make. 


availMoveSspice ■ AGENT x {AGENT POSITION) x {POSITION M) x {AGENT M) 
X {POSITION N) X {AGENT N) x {AGENT seq BIT) x N 
-^FPOSITION 

\/x, agent : AGENT] positions : AGENT h-» POSITION] vision : M; 
sugar, spice : POSITION N; agentSpice, agentSugar : AGENT hh- N; 
culture : AGENT h-» seq BIT • 

availMoveSspice{cigent,positions, sugar, agentSugar, spice, agentSpice, culture, vision) = 

{I : POSITION] X : AGENT \ I G distance{l,positions{agent)) < vision 
A positions{x) = I ^ {agentSugar{x)+SigentSpice{'x) < 
agentSugar{agent) + agentSpice(ag) 

A tribe{culture{x)) ^ tribe{culture{agent))) 

A {{distance{positions{x),l) < vision) (3) 

A tribe{culture{x)) ^ tribe{culture{agent))) ^ 

agent Sugar {x)+agentSpice{x.) < agentSugar{agent)+SLgentSpice{agent) 

Prewar d{l, sugar, positions, agentSugar, GO MB AT LI MIT)) 

+reward(l, spice, positions, agentSpice, SPICECOMBATLIMIT) • /} 
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— C OTTlbcit - 

ASugarScape 

step' = step 

maxSugar' = maxSugar 

maxSpice' = maxSpice 

pollution' = pollution 

loanBook' = population' < loanBook > {population' < {runloanBook)) 

spiceLoanBook' = 

population' < spiceLoanBook > (population' < (ranspiceLoauBook)) 

Wag : AGENT] I : POSITION • 

ag G population' ^ (3) 

{sex'{ag) = sex{ag) 

A vision'{ag) = vision{ag) 

A age'{ag) = age{ag) 

A maxAge'{ag) = maxAge{ag) 

A children'{ag) = children{ag) 

A agentCulture'{ag) = agentCulture{ag) 

A agentimmunity'{ag) = agentlmmunity{ag) 

A metabolism'{ag) = metabolism{ag) 

A initial Sugar' {ag) = initial Sugar{ag) 

A initialSpice'(ag) = initialSpice(ag) 

A spiceMetabolism'(ag) = spiceMetabolism(ag) 

A diseases' {ag) = diseases{ag)) 

{population', position', sugar', agentSugar', agentSpice') = 

applyAllCombatspice{asSeq{population),population,position, sugar, 
agentSugar, agentSpice, vision, agentCulture) 
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applyAllCombatspice : seq AGENT x F AGENT x {AGENT POSITION) 

X {POSITION M) X {AGENT >+^N)x {AGENT M) 

X {AGENT ^N)x {AGENT ^ seq BIT) 

-)■ 

{F AGENT X {AGENT POSITION) x {POSITION N) 

X {AGENT N) X {AGENT N)) 

Mhead : AGENT] tail : seq AG ENT; pop : F AGENT; 
positions : AGENT POSITION; sugar : POSITION N; 
agSugar, agSpice : AGENT h-s- N; 
vision : AGENT -H N; culture : AGENT -H seqBIT • 

apply AUG ombat{{), pop, positions, sugar, ag Sugar, ag Spice, vision, culture) = 
{pop, positions, sugar, ag Sugar, ag Spice) 


apply AUG ombat{{head) ^ tail, pop, positions, sugar, ag Sugar, ag Spice, vision, culture) = 
if {head G pop)then 

apply AUG ombatspice{tail, 

singleEightspice{head, pop, positions, sugar, agSugar, agSpice, vision, culture)) 

else 

apply AUG ombatspice{tail, pop, positions, sugar, agSugar, agSpice, vision, culture) 
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singleFight spice ■ AGENT x ¥ AGENT x {AGENT POSITION) x {POSITION >+^ 
X {AGENT N) X {AGENT N) x {AGENT ^N)x {AGENT ^ seqBIT) 

{¥ AGENT X {AGENT POSITION) x {POSITION M) 

X {AGENT N) X {AGENT N)) 

Vac/ : AGENT] population : F AGENT; positions : AGENT POSITION; 
sugar : POSITION N; agSugar, agSpice : AGENT h-> N; 
vision : AGENT -H- culture : AGENT seqBIT • 

singleFightspice{0'g, population, positions, sugar, agSugar, agSpice, vision, culture) = 
if {availMoves{ag,positions, sugar, agSugar, agSpice, culture, vision{ag)) = 0) then 
{population, positions, sugar, agSugar, agSpice) 

else 

3/oc : POSITION; available : F POSITION; 

'SotherLoc ; POSITION \ 
loc, otherLoc G 

availMoveSspice{ag, positions, sugar, agSugar, agSpice, culture, vision{ag)) 

A other Location ^ location 

reward{loc, sugar, position, agSugar, OOM BAT LIMIT) 

+reward(loc, spice, position, agSpice, SPICECOMBATLIMIT) 

> rewar d{other Loc, sugar, position, agSugar, GO MB AT LI MIT) 

+reward(otherLoc, spice, position, agSpice, SPICECOMBATLIMIT) 

{distance{position{ag), loc) > distance{position{otherLoc), position'{ag)) ^ 
reward{loc, sugar, position, agSugar, GOMBATLIMIT) 

+reward(loc, sugar, position, agSpice, SPICECOMBATLIMIT) 

> reward{otherLoc, sugar, position, agSugar, GOMBATLIMIT) 

+reward(otherLoc, spice, position, agSpice, SPICECOMBATLIMIT)) 

• 

{population \ {positions {loc)}, {positions ^ {loc}) © {ag i-A loc}, sugar © {loc i-A 0}, 
agSugar © {ag i—)■ agSugar{ag) + 

reward{position'{ag), sugar,position, agSugar, GOMBATLIMIT)}, 

agSpice © {ag i-A agSpice(ag) + 

reward(position'(ag), spice, position, agSpice, SPICECOMBATLIMIT)}) 
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_ SynchronousCombat spice _ 

ASugarScape 

step' = step 

maxSugar' = maxSugar 

maxSpice' = maxSpice 

pollution' = pollution 

loanBook' = population' < loanBook > {population' < {renaloanBook)) 

spiceLoanBook' = 

population' < spiceLoanBook > (population' < (ran spiceLoanBook)) 

population' C population 

sugar' = sugar © {p ; POSITION \ p G renaposition' • p i—)■ 0} 

Vap : AGENT] I : POSITION • 
ag G population' ^ 

{sex'{ag) = sex{ag) 

A vision'{ag) = vision{ag) 

A age'{ag) = age{ag) 

A maxAge'{ag) = maxAge{ag) 

A children'{ag) = children{ag) 

A agentCulture'{ag) = agentCulture{ag) 

A agentimmunity' {ag) = agenti immunity {ag) 

A metabolism'{ag) = metabolism{ag) 

A initial Sugar' {ag) = initialSugar{ag) 

A initialSpice'(ag) = initialSpice(ag) 

A spiceMetabolisni(ag) = spiceMetabolisni(ag) 

A diseases' {ag) = diseases{ag) 

A agent Sugar' {ag) = agent Sugar {ag) 

prewar d{position' {ag), sugar, position, agentSugar, COMB ATLIMIT) 

A agentSpice'{ag) = agentSpice{ag) 

+reward(position'(ag), spice, position, agentSpice, 
SPICECOMBATLIMIT) 

A position'{ag) G 

availMovesspice{ag, position, sugar, agentSugar, 
agentSpice, agentCulture, vision{ag))) 

) 

ag G population \ population' 

3x : AGENT • position'{x) = position{ag) A tribe{culture{x)) ^ tribe{culture{ag)) 

{I G availMovesspice{o-g, position, sugar, agentSugar, agentSpice, agentCulture, vision{ag)) 
A reward{l, sugar, position, agentSugar, COMBATLIMIT) 

+reward(l, spice, position, agentSpice, SPICECOMBATLIMIT) 

> rewar d{position' {ag), sugar, position, agentSugar, COM BATLIMIT) 

+reward(position'(ag), spice, position, agentSpice, 
SPICECOMBATLIMIT) 

A distance{position{ag), 1) < distance{position{ag),position'{ag))) 

=P 3x : AGENT • position'""{1) = x A position{x) ^ I 
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5.17. Credit. Credit is defined with one resource. It is incorporated into a dual 
resource simulation but no extra information is given as to what changes were 
made, if any. The most logical approach is to assume that spice loans are admin¬ 
istered the the exact same way as sugar. We create a separate system of loans for 
spice that is dealt with in the exact same manner as sugar loans. Our specihcation 
is now split into two parts, the parts dealing with sugar, already specihed, and the 
parts dealing with spice, which are copies of their counterparts. The amount of a 
resource available to be borrowed now depends on which resource we are talking 
about. 


amtAvailnew '■ N x SEX x N x M N 

V age, resource, baseAmt ; M • 
amtAvailnew{cige, male, resource, baseAmt) = 

resource div 2 age > MALEFERTILITYEND 
amtAvailnew{,CLge, male, resource, baseAmt) = resource — baseAmt 
{age < MALEFERTILITYEND A 
isFertile{age,male) A resource > baseAmt) 
amtAvailnew{o.ge, male, resource,base Amt) = 0 
{age < MALEFERTILITYEND A 
-1 {isFertile{age,male) V resource > baseAmt)) 
amt Avail new {o-g^-, female, resource, baseAmt) = 

resource div 2 age > FEMALEFERTILITYEND 
amtAvailnew{(^gG, female, resource, baseAmt) = resource — baseAmt 
{age < FEMALEFERTILITYEND A 
isFertile{age, female) A resource > baseAmt) 
amt Availnew{age, female, resource, baseAmt) = 0 
{age < FEMALEFERTILITYEND 
A -1 {isFertile{age, female) V resource > baseAmt)) 
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_ Pay Sugar Loans _ 

AAgents 

EStep 

population' = population 

sex' = sex 

position' = position 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentCulture' = agentCulture 

agentimmunity' = agent Immunity 

children' = children 

diseases' = diseases 

metabolism' = metabolism 

initialSugar' = initialSugar 

initialSpice' = initialSpice 

spiceMetabolism' = spiceMetabolism 

spiceLoanBook' = spiceLoanBook 

agentSpice' = agentSpice 

3dueLoans,newLoans : {AGENT {AGENT x (N x M))) • 

dueLoans = loanBook \> {Tan{loanBook) > {a : (M x N) | second{a) = step}) 

{loanBook', agentSugar') = 

payExclusiveLoans{chooseGonflictEreeSets{dueLoans), agentSugar, loanBook) 






104 


JOSEPH KEHOE 


_ PaySpiceLoans _ 

AAgents 

EStep 

population' = population 

sex' = sex 

position' = position 

vision' = vision 

age' = age 

maxAge' = maxAge 

agentCulture' = agentCulture 

agent Immunity' = agent Immunity 

children' = children 

diseases' = diseases 

metabolism' = metabolism 

initial Sugar' = initial Sugar 

initialSpice' = initialSpice 

spiceMetabolism' = spiceMetabolism 

loanBook' = loanBook 

agentSugar' = agentSugar 

3dueLoans,newLoans : [AGENT •H- [AGENT x (M x M))) • 

dueLoans = loanBook > [Tan[spiceLoanBook) > {a : (M x M) | second[a) = step}) 

[spiceLoanBook', agentSpice') = 

payExclusiveLoans[chooseGonflictEreeSets[dueLoans), agentSpice, spiceLoanBook) 






_ MakeLoans Sugar 

A Agents 
EStep 
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population' = population A sex' = sex 
position' = position A vision' = vision 
age' = age A maxAge' = maxAge 

agentCulture' = agentCulture A agentimmunity' = agent Immunity 
diseases' = diseases A children' = children 
metabolism' = metabolism 
initialSugar' = initialSugar 

initialSpice' = initialSpice 
spiceMetabolism' = spiceMetabolism 
agentSpice' = agentSpice 
spiceLoanBook' = spiceLoanBook 

3neujLoans : ¥{AGENT x (AGENT x (M x N))); 

V ag, lender, borrower : AGENT ; amt, due : M • 
loanBook' = loanBook U newLoans 
ag G dovanewLoans ^ 

agentSugar'(ag) = agentSugar(ag) — totalLoaned(ag, new Loans) 
ag G dom(ran newLoans) 

agentSugar' (ag) = agentSugar(ag) + totalOwed(ag, new Loans) 
ag ^ dom(newLoans) U dom(ran newLoans) 
agentSugar'(ag) = agentSugar(ag) 
willBorrow(age(ag), sex(ag), agentSugar'(ag), 

iaii(loanBook' n {a ; AGENT x (AGENT x (N x N)) 

I borrower (a) = borrower (loan)})) 

-1 3ag2 : AGENT • canLend(age(ag2), sex(ag2), agentSugar'(ag2)) 

A adjacent(position(ag2),position(ag)) 
total Loaned(ag, new Loans) < 

amtAvailnew((^9^(0‘9)y sex(ag), agentSugar (ag), GHILDAMT) 
totalOwed(ag, new Loans) < GHILDAMT — agentSugar (ag) 

(lender, (borrower, (amt, due))) G newLoans 

(canLend(age(lender), sex(lender), agentSugar(lender)) 

A willBorrow(age(borrower), sex(borrower), agentSugar (borrower), 
{borrower} < (rasi loanBook)) 

A amt < 

min({amtAvailnew(o-ge(lender), sex(lender), agentSugar (lender), 

GHILDAMT), GHILDAMT — agent Sugar (borrower 
A due = step + DURATION 
A adjacent(position(lender),position(borrower))) 
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_ M akeSpiceLoans _ 

ESpiceLattice 

ASpiceAgents 

EStep 

population' = population A sex' = sex 
position' = position A vision' = vision 
age' = age A maxAge' = maxAge 

agentCulture' = agentCulture A agentimmunity' = agent Immunity 

diseases' = diseases A children' = children 

metabolism' = metabolism 

spiceM etabolism' = spiceMetabolism 

agentSugar' = agentSugar 

loanBook' = loanBook 

initialSugar' = initialSugar 

initialSpice' = initialSpice 

3 new Loans : F{AGENT x {AGENT x (M x N))); 

V ag^ lender, borrower : AGENT ; amt, due : M • 
spiceLoanBook' = spiceLoanBook U new Loans 
ag G dom new Loans ^ 

agentSpice'{ag) = agentSpice{ag) — totalLoaned{ag, newLoans) 
ag G dom(ran newLoans) 

agentSpice'{ag) = agentSpice{ag) + totalOwed{ag, new Loans) 
ag ^ dom{newLoans) U dom(ran newLoans) 
agentSpice' {ag) = agentSpice{ag) 
willBorrow{age{ag), sex{ag), agentSpice'{ag), 

iaxi{loanBook' n {a ; AGENT x {AGENT x (N x N)) 

I borrower {a) = borrower {loan)})) 

-1 3ag2 : AGENT • canLend{age{ag2), sex{ag2), agentSpice'{ag2)) 
A adjacent{position{ag2),position{ag)) 
total Loaned{ag, new Loans) < 

amtAvailnew{o-ge{ag), sex{ag), agentSpice{ag), SPIGEGHILDAMT) 
totalOwed{ag, new Loans) < SPIGEGHILDAMT — agentSpice{ag) 

{lender, {borrower, {amt, due))) G newLoans ^ 

{canLend{age{lender), sex{lender), agentSpice{lender)) 

A willBorrow{age{borrower), sex{borrower), agentSpice{borrower), 
{borrower} < {ran spiceLoanBook)) 

A amt < min{{amtAvailnew{(^ge{lender), sex{lender), 

agentSpice{lender), SPIGEGHILDAMT) 
SPIGEGHILDAMT — agentSpice{bor rower))}) 

A due = step + DURATION 
A adjacent{position{lender),position{borrower))) 
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5.18. Rule Application Sequence. 

T tck spice 

[9 Growbackspice I 9 SeasonalGrowbackspice] 

[9 Ad OVCimCTltijCLsicSpice I 9 pollutions pice 9 d^olltlttOTlJDt f f tlStOTispice) I9 ^ spice] 

{9 Inheritance spice} is Deathspicels R&plctcement spice I9 1^9 ent Mating spice]} 

{9 Gulture}}^ Pay Sugar Loans 9 PaySpiceLoans 9 MakeSugar Loans 9 M akeSpiceLoans} 
{9 Transmission spice 9 ImmunePesponsCspice} is Trade} 


6. Conclusions 

We have shown that it is possible to apply formal methods fruitfully in the 
area of ABSS and produced a full formal specihcation of the Sugarscape family of 
simulations. It is, to the best of our knowledge, the hrst formal specihcation of the 
entire Sugarscape simulation family. The purpose of the specihcation is to provide 
a clear, unambiguous and precise dehnition of Sugarscape. The specihcation has 
identihed many ambiguities and/or missing bits of information in the original rule 
dehnitions. Where there is an obvious way of removing these ambiguities we have 
done so. If there is more than one possible solution we have identihed them and 
chosen the most likely one. 

The issues with the model dehnition that we have encountered can broadly be 
grouped into three main types: Lack of Clarity, Missing Information and Sequential 
biases. 

Lack of Clarity: The rules, although simply stated in the appendix, lack 
clarity in their dehnition. Only one version of each rule is presented even 
when many variations are referred to in the text. The variations presented 
cannot always be used together, for example the Movement rule dehned in 
the appendix is not the variant required if the pollution rule is also used. 
Our specihcation brings them all together in one place for ease of reference. 

Missing Information: Missing or incomplete information is the biggest 
cause for concern. In many cases we can work out the most likely answer 
based on context but in some cases there is not one dehnitive correct an¬ 
swer. If there was more than one arguably correct solution we chose the 
simplest. How we hll in these blanks can have a big ehect on how the 
simulation proceeds. These ehects may be important if we are trying to 
compare diherent implementations of Sugarscape. If we take the disease 
transmission rule, for example, questions that are unanswered include: 

(1) Once an agent gains immunity from a particular disease, do we re¬ 
move that particular disease from the set of diseases that the agent is 
carrying, or is the agent still a carrier? 

(2) When we transmit a disease do we only transmit diseases that we carry 
and have no immunity for, or, can any disease we carry be transmitted? 
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(3) The Mating rule omits important information about parents contribut¬ 
ing half of their resources to their offspring. This has a huge effect on 
how mating works in a simulation. 

By replacing each ambiguous interpretation with one simple and precise 
interpretation we allow different developers to replicate their results and 
benchmark them against each other. All hidden assumptions that could 
serve to advantage one implementation over another are excised. 

Sequential Biases: Sugarscape is based on the assumption that it will be 
implemented sequentially. While this may have been a good assumption at 
the time it was written it is not now necessarily the case. Improvements in 
processing speed have recently been attained mainly through the introduc¬ 
tion of concurrency. Simulations are now almost always run on multicore 
or even multiprocessor machines. 

The Z specihcation is free from assumptions about implementation. It 
achieves this without having to specify or constrain in any way what conflict 
resolution or avoidance strategies are employed. This leaves developers 
the freedom to try out different approaches as suits their implementation 
platform. 

Because the specihcation is high level and only dehnes the before and after 
state of each rule it makes few assumptions as to how any rule will be imple¬ 
mented. All inherent biases towards a sequential implementation are removed. 
Implementers have complete freedom as to what programming model they em¬ 
ploy (Object-oriented, imperative, functional, or any concurrent approach). Any 
simulations that adhere to the standard can be properly compared in terms of 
performance or patterns of behaviour. This will put on a hrmer foundation any 
claims made by researchers about their implementations. 

6.1. Further Work. Further work remains to be done in getting agreement from 
the ABM community on the decisions made in producing this interpretation of 
Sugarscape. Any incorrect assumptions made in removing ambiguities need to 
be identihed and agreed upon. This provides a route to address the issues of 
replication of experimental results in ABSS. 

Sugarscape can now be used as a benchmark (or rather set of benchmarks) for 
ABM implementers. This is particularly useful for those proposing new approaches 
to concurrency that promise performance improvements. Current trends, for exam¬ 
ple, include the use of Graphics Processor Units (GPUs)[Deissenberg et ah, 2008, 
Lysenko and D’Souza, 2008, Richmond et ah, 2009], containing hundreds to thou¬ 
sands of individual processors. These approaches tend not to be tested on the 
more complex rules in Sugarscape (such as Combat, Inheritance and Trade) as 
they are not easily parallelized. By providing a precise and full set of these rules 
it is now possible for researchers to properly compare how different models cope 
with more complex and more realistic ABMs. 
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Z itself is rather verbose and can be hard to parse when reading. The lack 
of modularity made function dehnition signatures overly long. The available tools 
such as CZT [Malik and fitting, 2005] make the process of writing the specihcation 
easier but I have altered the specihcations to remove bracketing where I thought 
it made the specihcation easier to read even if this was hagged as an error in the 
type checker. These issues could be overcome through the use of a variant of Z 
such as Object-Z or Alloy. The issue of whether ABM modellers would be willing 
to use formal specihcations remains unknown. 

There are diherences between the outcomes of the synchronous and asynchro¬ 
nous approaches. Sugarscape assumes an asynchronous approach and this ahects 
the style of specihcation that we use. We have shown in the case of combat the 
diherences in a synchronous and asynchronous specihcation. While we regard the 
synchronous specihcation as somewhat simpler to produce but others may disagree. 
We tackle the question as to which approach is the more correct elsewhere. 
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